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computes infrared radiances for an earth's limb viewing geometry. fhe 

nominal spectral region of this study lies between 2.7 and 25 "1.ro- 

meters and emphasis is placed on radiation originating at alt.tudes 

between 70 and 500 kin. The physical model is described, with emphasis J 
on the changes required in extending its usefulness. Application of the 

computer program is described. Model predictions are compr-ed with radi- 

ance data obtained in the ICECAP auroval studies program. inis compari- 

son is used to illustrate uncertainties in results due to ascumptions 
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1.0 INTRODUCTION 


The research reported here has had as its major objective the develop- 
ment and improvement of infrared radiance and transmission models of the 
atmosphere. These models are designed to simulate the radiative properties 
of the atmosphere to provide prediccions for Air Force and other Department 
of Defense design and surveillance programs. Comparison with available ex- 
oerimental measurements provides greater understanding ot the atmosphere and 
serves to check the adequacy of the models. 

The main areas of effort have been 

1. Extension of the applicability of a current 
model on non-equilibrium atmospheric infrared 
radiance between 2.5 and 25 inicrons. 


2. Development and improvement of thermal equi- 
librium mod*ls of the infrared transmission 
and radiance for the earth's lower atinosphere. 


The goal of the first area has been the further development of a com- 
puter program to simulate the natural infrared radiance background of the 
earth's upper atmosphere. The nominal spectral region under study lies be- 
tween 2.7 and 25 micrometers and emphasis is placed on radiation originating 
at altitudes between 71) and 500 kilometers. The general problem area is of 
interest for systems design, military surveillance and the advancement of 
knowledge about physical processes in the upper atmosphere. The immediate 
application of this work is to aid in developing optimum invrared background 
measurements programs and in interpreting the results of such measurements. 

This work is an extension of the study of Corbin, et al (1969) and 
Degges (1972, 1974}. The former investigated the natural infrared 
backaround of the earth in 5 to 25 micrometer spectral reqion, with 
the qoal of estimating earth timb viewing radiances for tangent heights from 
the surface to 500 km altitude. For convenience, their study divided the 
atmosphere into two regions with a division at 70 km. Below 70 km the atmo- 
sphere was assumed to be in thermal equilibrium. Adove 70 kn axplicit cal- 
culations were made cf processes which excite and de-excite muiecuia: vibra- 
tional and rotational levels which are the soucce of infrared radiation. 


a8 cf vale natal gO 6s 


Their study concentrated on radiation from water vapor, carbon dioxide, 
Ozone, Nitric oxide and nitrous oxide, which are principal radiating species 
in the spectral region considered. In addition, nitric acid was included in 
the lower atmosphere work and estimates were made of radiation to be expected 
from particulate matter suspended in the atmosphere. 

Corbin, et al (1969), presented models for the lower atmosphere for 
a wide range of seasonal and latitudinal conditions. This was not posssible 
for the abundances of most minor neutral species. More data has since become 
available, particulurly for nitric oxide and the hydroxy) radical, but at 
present it appears that the best means of estimating abundances of important 
infrared emitting species is chemical rate equation integrations including 
molecular diffusion and eddy mixing. Even calculations involving transport 
properties cannot always be accepted because the values of eddy mixing co- 
efficients sre to a large extent only informed quesses and published calcula- 
tions often use outdated rate coefficients. Degaes (1972) reported a comouter 
program with which to determine diurnal variations in abundances of minor 
species, to investigate the effects of changes in assumed eddy mixing cuef- 
ficients, and to estimate the effects of new determinations of chemical rate 
coefficients. 

A second area of study necessary to improve the radiance model includes 
the physical arocesses that control the population of infrared emitting states 
of atmospheric molecules. Except for the cure curctional radiation from mole- 
cules such as water, the degree of cacitation of vibrational levels determines 
the radiation from infrared emitting molecules. The most important mechanisms 
are collistonal exc‘ tacion and de-excitation and absorption and re-enission of 
electromagnetire radiation, 

In the troposphere and lower stratosphere, collisional processes are 
rapid enough to control the population, of vibrational levels. Above 30 to 
5d Am, Nowever, collisvonal eacitation becomes less efficient and radiative 
processes become important. The combined effects of collisional and radtative 
prevesses must therefore ve considered. Below an altitude of about 9 he 
molecular nitrogen and vaygen are the most important callisive pat iners, 

Above that altitude, atome oxygen becemes duper tant. beta in esciting attiay 


oxide ana in determining malecular vaygen and nidiegen widratie uel Temperatures, 


. 
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In determining the effects of radiation on the populations of vibra- 
tional levels of infrared active molecules, it is necessary to separate the 
radiation of a single change in vibrational quantum numbers from the rest of 
the radiation field. The previously reported studies did this by assuming 
a Doppler line shape for the individual rotational lines of a band and were 
abla to obtain aceqvate numerical approximations for radiative transfer func- 
tions appropriate to single bands of linear molecules, and less accurately, 
for water vapor and ozone bands lying in the spectral region of interest. 

The principal areas of research and program revision reported by 
Degges (1974) included a further study of band radiance modelling and the 
effects of the line shape on radiance computations. In addition to some 
exploratory calculations using a Voigt line profile instead of a Doppler 
line shape, some asymptotic analytical results were obtained which provide 
an estimate of the accuracy of the numerically cotained band functions at 
large optical thicknesses. It was found that when a Voigt profile is used, 
the connection between the band models and functions for a single jine fs 
more readily seen than when a Doppler line shape is used. This basic mate- 
rial on radiative transport is reviewed in Section 2, and additional results 
are presented, 

The status of current knowledge of atmospheric composition is re- 
viewed in Section 3, Recent measurements of chemical rate coefficients 


have resulted in changes in estimates :f abundances of some molecular species. 


An attempt is made to set upper and lower bounds for mixing ratios at high 
altitudes. 

Resiits of calculations are given in Section 4. We discuss the ef- 
fects of changes made in the program and compare model predictions with 
ICECAP radiance measurements. This experimental data permits us to better 
assess the validity and limitations of the model. 

In Section 5 we review the limitations of low altitude models which 
assume thermal equilibrium for vibrational populations and estimate the 
altitudes below which such models are valid. 

In Section 6 we discuss the high resolution Vine-by-line program, 
This is treated more fully in a separate scientific report. The computer 
programs are documented in Section 7. 


Mater at 


mano tare 


RE EE RO SE IIE RENE SS RR ne Ein ~ it Sil acces +” cP epee waitin hebenise eat Wewatich haan tie T= aeths Sapte ee ae ate paid mertietenende net SF yn a che! 


2.0 RADIATIVE PROCESSES AND MODELLING 
2.1 Introduction 


The physical processes and assumptions used in developing a 
model for transport of infrared radiation in the upper atmosphere are dis- 
cussed below, beginning with a general overview and discussion of a simpler 
problem. Next, the interaction of electromagnetic radiation with an isolated 
atom or molecule is described. An understanding of the basic processes in- 
volved is the first step in any computation of the radiation from a gas or 
mixture of gases. Subsequently, absorption from a continuum source, a spec- 
traily varying source and fluorescent processes are treated. The effects of 
optica? thickness when the path of the radiation includes a large number of 
molecules of the same species is considered for spectral lines with a Voigt 
profile, which includes Doppler (Gaussian) and Lorentz shapes as limiting 
cases. Finally, a model for the transport of radiation for a single vibra- 
tion-rotation band is described, and the effects of temperature estimated. 

Figure 1 illustrates the gecmetry employed in the current 
modelling study. A molecule at Point P absorbs radiation from three sources: 
(1) a collimated beam from the sun, (2) the lower boundary which radiates as 
a black body, and (3) other molecules betwecn the lower and upper boundaries. 
Radiative transfer functions are needed to compute the effects of absorption 
of radiation by other molecules between each of the sources and the molecule 
at Point P. 

In the absence of collisional excitation and de-excitation, 
the simplest integral equation that may be written to describe radiative trans- 
port ina finite plane parallel atmosphere is 


tT 
ata) = Ff alle = th) oleae + (0 (a) 
tY) 


This equation assumes that there is only one radiative transition involved 
and that the radiation emitted by a molecule after absorption is independent 
of the direction of incidence. Here t is a measure of the optical path length 
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Finure 1. Geometry assumed in der‘ving the radiative transfer functions. 
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in the vertical direction. J(t) is a dimensionless measure of the population 
of the upper, emitting molecular state. K,(t) is the radiative transfer func- 
tion which accounts for absorption of radiation emitted by molecules at a 
vertical path length t distant from the absorbing molecule. S(+r) is an 
auxiliary function which includes the effects of radiation originating outside 
the atmosphere. If there is only a blackbody source located at the lower 
boundary from which optical path length is measured (t = 0 at the lower 
boundary) then 


S(t) = Ny Ko(t) (2) 


Here No is the radiance of the lower boundary in arbitrary units such as 
watts cm“-ster”! or the spectral radiance in units such as watts em? cm™2. 
ster”! (wavelength interval)7?, 

In the gray or spectrally independent case, when Equation (2) 
holds, the kernel function Ky (1) is the first exponential integral Ey (1) and 
the function Ko(t) describing the absorption from the blackbody source is the 
second exponential integral E(t). These functions are described in detail 
in such works as Kourganoff (1963). A complete analytic solution {is available 
and the auxiliary functions required for a numerical solution have been tabu- 
lated by King (1956). Figure 2 plots the solution of Equation (1) when Equa- 
tion (2) holds for fie gray case for values of total optical thickness Tage 
0, 0.1, 1.0, anu ad. When the total optical thickness is small, the effective 
excitation of the molecule at any poi~t in the atmosphere is half that of a 
molecule at the (vibrational) temperature of the blackbody. As the total 
optical thickness increases, the excitation of the molecules near the black- 
body inc: cases, while that of molecules near the top of the atmosphere de- 
creases. This occurs because as the otpical thickness increases, molecules 
near the lower boundary absorb an increasing amount of radiation reradiated 
from molecules above them. At the same time, less radiation penetrates to 
higher altivude. The physics of the problev. require that if a second black- 
body of the same temperature be placed at the top of the atmosphere, the 
population of radiators at all points would be uniformly that corresponding 
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Figure 2: Effect of optical thickness on relative vibrational 
excitation of atmosohere illuminated from below by 


black body source (Gray Case). 
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to the blackbody temperature. This symmetry and the linear character of the 
integral Fquation (1) results in the relation d(:) + Lee -t)- No: 

When a collimated beam of light is incident on the top of the 
atmosphere, another term must be added to the function S(i) in Equations (1) 
and (2) to account for this. Tf the flux is incident at a zenith angle with 
cosine yond has an intensity in arbitrary units such as nF watts ew"! normal 
to its direction, the term added to S(t) has the form FM, eter -~a)/u). 
In the aray case, Mr) is simply exp(-1). The corresponding functions for 
individual spectral Pines and bands will be described after a closer Took at 
the interaction af efectroamagnetic radiation with a single molecule. 


ocr Interaction of Radiation With A Single Molecule 


When treating @ probler in which the number of photons enitted 

or absorbed by an atom or molecule mist be considered rather than total 
eneray eitted, dt ts conventent to examine the physics of the interaction 
Of radiation with matter through use af the Einstein A and 8 coeffictents, 
In treating a radiative problem such as the current study of high altitude 
infrared radiances which alsa involves transfer of energy between molecules 
through collisdons, this appraach is a necessity. In the following discus- 
ston. caus units are used except where oxaplicitly stated, 

the binstein Ave coetficient ts the rate at which an atom or 
mwlacule ds an excited state spontaneously emits radfation. tt has the units 
photons ?secmotecule, Tn an assemblage of molecules all inoa state which can 
radiate. the average rate of radtation fs given by the Einstein Aue cet - 
ficient for the particular radfative transition from a state of higher: eneray 
(subscript ud Cooonme of lower eneray (subscript ¢). 

In addition toa spontaneous emission, Gea other processes accur, 
These are absorption, chara: tertced by the Einstein coetfictent for absorp- 
tion Reus and tudtced eatssdan, characterized by the Eiastein coeffictent 
tor induced amiss ton Bie Yyoa Blackbody radiation field with spectral 
density, OTE the probability fhat an atan or mofecule abserbs a quantum 
Of Madiattan on unit Time is Reha! When radiatien ad is expiessed 
COS UNETN, ers ses cay the coet fay done se has undgs ch Serge sec. Sim 


flardv. an atom or amelecule moa radiation Cfedd with the above spe (ral 
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density and in an excited state emits radiation at the rate BigP (%,) in ad- 
dition to that characterized by the coefficient Ae: The Einstein coef- 
ficients are properties of an atom or molecule. If one coefficient is known, 
the others may be determined by the relations 


3 
Ae * BnhB /A (3) 


and 
pb * Sub uy (4) 


Thus, it is possible to treat the problem knowing only the 
Einstein Aue coefficient. 

Experimental determinations of the interaction of infrared 
radiation with gases are usually made by measurements of light absorption. 
Experimental results may be expressed as integrated absorption from a con- 
tinuum source, 


stow N Aa wy oe eee ~ fe (S) 
ty 1 Br ac 


or as line or band strengths 


Say = S'gyéPe (6) 


where p is the pressure. 

In the visible and ultraviolet regions of the spectrum, it 1s 
customary to use oscillator strengths to characterize the interaction of 
radiation and matter. In this case, the defining equation is 


The resulis of theoretical computations of the interaction cf 


matter with radiation are frequently expressed as dipole moments, Moe: These 
are related to the Einstein A coefficient by 
2 
4g 
641 g 
Ay? ya a [Mh (8) 
uk 3hr 9, uk 


The rate at which a sinale atom or molecule absorbs and reradiates light from 
a parallel beam from a continuum source at a wavelength do corresponding to 
a resonance transition is (Mitchell and Zemansky. 1961): 


2 
» Me 1 ¢ 42 
G : Lofors (photons/sec) 
(9) 
. =13 2 
8.852 x 10°" 1f,A5  (photons/sec) 


where ly is the continuum source flux and ft is the oscillator strength. 

If the incident flux is expressed in photons /cm*-ym, and oscillator strength 
is replaced by the equivalent line or band strength S expressed as cm™“atm7?, 
then 


@= 3.72 x 107° 1 sa? (10) 
Implicit in this expression is the convenient relation between the Einstein 


Ae coefficient and line or band strength with wavelength dimensions ym 


Aig * 2804S (om™atm™*)/a? (um) (un) 

When the source is an extended blackbody surface with a spec- 
*ral radiance Ny» the G factor is obtained by integration over the solid angle 
subtended by the source. For radiation from a plane parallel surtace of in- 
finite extent, 


23 sa (a2) 


G = 1.169 x 10° 

Another purely radiative mechanism of importance in upper 
atmosphere radiance computations is fluorescence. A striking example occurs 
for the weak 10.4 um C0, band. Absorption of solar radiation by C0, 4.3 um 
band maintains the V3 vibrational mode at a vibrational temperature of about 
290 K. During the day, this leads to an enhancement of emission in the 
10.4 um band by as much as two orders of magnitude when viewed along the 
earth's limb from outside the atmosphere. 

Fluorescence occurs when the upper state of the resonance 
transition reradiates into other than the initial states. The relative 
probability of a radiative transition to a lower state, n from the upper 
state, u when there are k possible lower state is: 


Aun 
at 
yw 


P 2 


un (13) 


where the A, are the Einstein A coefficients for the spontaneous emission 
probabilities for a transition from the upper state to the jth lower state. 

In the case of molecular resonance-fluorescent scattering at 
low temperature, only the lowest vibrational level (v" = 0) in the ground 
electronic state is populated significantly. Absorption of a photon excites 
the molecule to a vibrational level, v' in the same or a higher electronic 
State. Radiation is possible from this siite to the initial level (v" # 0) 
or other vibrational levels (v" = 0) in the ground electronic state. The 
rate of resonance fluorescent scattering in a given electronic vibrational 
transition may be obtained by combining Equations (9) and (13) and intro- 
ducing appropriate subscripts: 


G,.. = 8.853 x 1077 1, £, 2 Avy (14) 
viy" vio vio v'o } Av'y" 
y" 
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Here the subscript v'v" refers to a transition from an upper state with vi- 
brational level v' to a lower state with vibrational level v". The subscript 
v'o refers to absorption of a photon by a molecule in the lowest vibrational 
level of the ground state, exciting the molecule to a vibrational leves v' 

in the same or higher electronic state. 


2.3 Radiative Transfer Functions for Lines and Bands 


In developing a realistic model of the interaction or radiation 
with matter, it is necessary to take into account the details of absorption. 
Line shape must be considered because it determines the amount of radiation 
escaping from an optically thick layer of gas. Thus, for large optical thick- 
ness (greater than a magnitude of about 10), the total emission from an iso- 
lated spectral line with Lorentz shape is proportional to the square root of 
the optical thickness, while that from a line with Doppler shape is propor- 
tional to the square root of *he logarithm of the optical thickness. In 
addition, absorption of radiation emitted at one place in the gas by a mole- 
cule at another location cannot be computed without knowledge of the line 
shape. 

At high altitudes, the profile of the spectral absorption curve 
of an individual line in a vibration-rotation band is predominantly due to 
velocity or Doppler broadening. Kuhn and London (1969), who investigated 
the heat budget of the atmosphere between 30 and 110 km, report that radiative 
transfer calculations require the use of the Voigt profile (including both 
collisional and Doppler broadening) at altitudes between 20 and 70 km. Altern- 
atively, the Doppler profile is adequate at higher altitudes where the Doppler 
shape is the largest contributor to line broadening. The original computational 
model therefore assumed a Doppler profile. 

The discussion below assumes a Voigt profile for a single line, 
and it will be shown that for small optical thicknesses and low pressures, 
the difference between calculations with Doppler and Voigt profiles is small. 
Asymptotic values for large optical thicknesses are quite different, however. 

The radiative transfer functions required to treat transfer of 
resonance radiation in a single line are well known for the Voigt profile 
and its two limiting cases, the Doppler and Lorentz profiles. Avrett and 
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and Hummer (1965) ‘eview the properties of four of these functions, giving 
series expanstons 1 or small values of optical thickness and asymptotic expres- 
sions for large values of the optical thickness. The notation used for these 
functions varies widely. That used here is adapted from that used by Ivanov 
and Shcherbakov (1965 a,b) in their tabulation of functions for the Doppler 
profile. 


The best known of these functions is probably that relating 
to the equivalent width of a line in absorption or to the total radiation 
emitted in a line from a source in thermal equilibrium at a constant temper- 
ature, the integrated absorptance, defined as 


L(a,t) = =) ° : on *flasx)) dx (15) 


Here f(a,x) is used to denote the Voigt profile, 


é 


re -t 
f(a,x) = af ra See (16) 


The parameter a is proportional to the ratio of the Lorentz half width at 


half maximum db and the Doppler half width at half maximum bp» 


b 
a= gt (1n2)!/? (17) 
1) 
The parameter x is proportional to the difference between the spectral fre- 


quency of the line center Vo and that at any other point in the line profile, 
v 


2(v - v.) 
x= — (an2)!/2 (18) 


The quantity 1+ is the optical thickness at the center of a line with Doppler 
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broadening alone, 


2 
éne” Nf \/2 
5 Sean (tn2) (19) 
mec by 


where N is the number density of ground state molecules and & {s the path 
length. The actual optica? thickness at the center of a line with Voigt 
profile is the product of the + of Equation (19) and exp (a2) times the 
co-error function of the parameter a, erfc(a). In general, Equation (15) 
must be evaluated numerically. For small values of +, a series expansion 
may be obtained in the form 


L(a, t) = - y {er 9%) (20) 
n=) 
where 
on(0) = 2 f fran]! (21) 


For all values of a, o,(a} is one. When @, (a) is divided by erfc a, a smooth 
variation between a = o (Doppler profile) and a = © (Lorentz profile} is ob- 
tained as shown in Figure 3 for small n. For large values of +, asymptotic 
series are available, as given for example by van Trigt (1968). 


1/2 
L(o,t) = 2 (tt) E + 2088 $53 ] (22a) 
1/2 2 
L(a#o,t) om /4 [ + a aa oof ] (22b) 


Equation (10) gives the rate at which a single molecule absorbs 
radiation from a parallel beam. If the molecule is located within an optically 
thick medium, the rate of absorption of radiation is reduced by a factor which 
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depends on the optical thickness, 


M (ast) = uf F(ayx)en tf lax) dx (23) 


For small values of t, a series expansion is 


© n 
M, (a,t) * ato fx) 944 (2) (24) 


Asymptotic expressions for large values of r are 


| 0. 2886 
Nits) eee Ws ee (25a) 
1 TET ant 
where 
M,(az0,t) 2 ey H/2 |, fi a- 2a? , (25b) 
1 al w/a T l oS at ears 


The function My is the derivative with respect to 1 of the function L. In 
addition to its use in modifying the G rfactor of Equation (10), it gives the 
probability of escape of a photon in a direction with optical thickness +t. 

A third function is defined as 


My(a.t) = al Fasx)® en TF(BX) gy (26) 


and has a series expansion 


%0 n 
Hare) =D fer) g ta) (27) 
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This function is proportional to the rate at which a molecule absorbs radi- 
ation from another molecule at a optical distance 1. This function would be 
used as the kerne] function Ky in Equation (1) if one dimensional radiative 
transport were being treated. It is included here because of its use in gen- 
erating the kernel function fer a plane parallel geometry. 

When the source of radfation is an infinitely extended plane 
blackbody, the function which is used to determine the rate at which a mole- 
cule absorbs radiation is given by 


Nyolast) = af f(a.x) E, [rF(a.x)] dx (28a) 
7 dt 

sof mie) $ (280) 
t t 


s/2 
* [ a (is) sinede (28c) 


The series expansion for this function is given by 


No last) ® $, (a) + (yt - t+ rent) e(2) (29) 
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Here a new integral is defined, 


2 
@ (a) = +f [f(a.x)} tn flax) dx 


and y is the Euler-Mascheroni constant 0.5772... The asymptotic expansions 
Ny 2 ( 0 ot : = Ties i ese (0 . 5 onl 0093 dare | ( 3 0a ) 
var “tnt “ } 
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may be obtained from Equations (29b) and @5). 
The kernel function for plane parallel geometry is given by 


a 2 
No (a.t) = ae [f(a.x)]} E, i fla.x)] dx (31a) 

¥T 

f M,(t) “ (31b) 
T 
Me 

F i stnede te} 

; Mo cos! case (23e) 


The series expansion useful for small 1 is 


oo n 
Noptaat) = 2 4 tnt) o, - o - y (~1). t (32) 


qe bent One 


and the asymptotic expansions are 


: 0193 
No (Ost) 2-5 ep i (0.5 eg Pek ‘) (33a) 


The recurrence relation 
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carries over into the corresponding band functions and provides a partial 
check on numerical calculations. 

In extending this treatment of radiative transfer to a vibra- 
tion-rotation band, in addition to assuming that each line hus the same Voigt 
profile shape, the following assumptions are made: 


1. Rotational levels witn.n a vibrational level 
remain in equilibrium with the translational 
temperature. This is valid up to at least 
150 km according to Goody (1964). 


“3 


there is complete redistribution of frequency 
within a single line following absorption of 
radiation. 


3. There is complete redistribution of aosorbed 
energy among all rotational levels. 


4. Line strenoths are used as if each line were 
at the band center, neglecting wavelength 
variation across a band, 


5. There {s no overlap of adjacent lines. 


6. Temperature variations at different levels do 
not affect the radiative transfer. 
Let 9) be the fraction of total band strength associated with a 
single line so that 


fy * 1 (35) 


and the sum over J may involve complete sunmation over rotational quantum num- 
ber for two or more branches. The optical thickness at the center of a single 
line will be given by TG) where + is the optical thickness that would result 
if the entire band strength were in a single line. Then the radiative trans- 
fer functions for a band may be defined as functions of radiative transfer 
functions for a single line. First, 
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rast) « J Ulas x99) (36) 


where the superscript b refers to an entire band, 
In orde™ that the correct derivative relations 


Nast) ge Care) (37a) 
yo (a.t) w= £ M>(a,t) (37b) 


and the integral relations of Equations (28b) and (31b) hold, 


wast) © j 9) Myles 199) (38a) 
nO(a.t) » 9% Malas t9y) (23b) 
NPotast) . j Gy Nyolas v9) (38c) 
M. (ast) . ’ 9 Noy (as 195) (38d) 


White these functions must be evaluated numerically, it is pos- 
sible to perform this evaluation in two different ways for small values of 1. 
Making the definitions 


G,° \ 9) (39) 
ee 
G ; 9) tn gy (39b) 
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The following series expansions are obtained: 
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NP (ast) = Go, (a) + Gye, (a) rent + Exo 
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For large values of the optical thickness parameter t, there 
is a sharp distinction between the asymptotic expansions for Doppler line 
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shape. When Doppler line shape is assumed, the asymptotic forms for the 
functions discussed above differ from those for a single line by a factor 
proportional to (anr)t/2 in the case of bands for a linear molecule. The 
numberical results suggest that the asymptotic expansions for an asymmetric 
rotor are multiplied by an additional factor proportional to tnz. In the 
case of a Voigt profile, there is no such change in the asymptotic form. 
The L functions retain their We dependence, Nie functions retain the + 
dependence, and the Mo and Noy functions retain the wile dependence. The 
chief complication with the functions for Voigt profiles is that for inter- 
mediate values of t the behavior is similar to that for a Doppler profile, 
while at some value of +t dependent on the parameter a, the functions begin 
to tend toward the asymptotic form. Similar behavior has been found by 
Malkmus (1967, 1968) for curves of growth for various random band models. 

A strict derivation of the asymptotic expansions for the 
Doppler profile functions has not been obtained. In the high rotational 
temperature limit, the first term in the expansions for parallel and per- 
pendicular bands of linear molecules appear to be 


Parallel Perpendicular 
(0,1) ~ vie Qp ent ; “ne Qp gm (41a) 
MP(o.t) ~ vie Qp/t 3 VRE Opt (41b) 
Mo(.t) ae Q,/x? 3 Vic Qp/t® (41c) 
L(o ‘ie le I, ent 3 in Q, tnt (41d) 
ae 2 R gr AR 
nO (o,t) - 1 Yae Q,/t 3 Vie Qo/1 (4le) 
yee eae) R 4 R 
ND (our) . 7 vue Q,/t? i Vie Qp/? (41f) 
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Here Qp is the rotational partition function and « is the he® /kT, the dimen- 
sionless ratio of rotational) quantum energy to kinetic temperature. 

The first term in the asymptotic expansions for Voigt profile 
band functions is easily obtained, as is illustrated for the L>(a,r) function. 
At large optical paths, L(a,t) for a single lina is from Equation (21b) ap- 
proximately ant/4 (az) t/2, Since the optical thickness at the center of a 
band is simply Gyt then 


carr) an7t/4 gl/2 (x9, )8/2 


or, denoting j 9° by 61/2 
b 2 1/2 
u'(a,r) ~ (at)“/* & (42a) 
nia 1/2 
and 


V2 

Wla.r) ~ La 81/2 (42b) 
V2 

Nolet) ~ ara (42e) 


w t 
Lar). bat ahi 872 (42d) 
3 3n / 
1/2 
b 2 a (42e) 
n° (a,t) - gil? G (42f) 
aver) ~ Sa7e37e Sie 
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In order to simplify programming, a single set of functions 


is used for all parallel bends of linear molecules and a single set is used 
for perpendicular bands. These functions vary with the rotational tempera- 
ture and in the program, in order to partially account for temperature vari- 
ation, the parameter 1 is not used, but rather a derived parameter To)? which 
is the optical thickness at the center of the strongest line in the band. 
This is either an R branch line for a parallel band or a Q branch line for 
a perpendicular band. 

In practice, rather than formally integrate a modified 
Equation (1) including a t dependent A to account for collisicnal de-excitation, 
and an additional source term to account for collisional excitation, the com- 
puter program sets an initial vibrational population at each alritude and 
computes the total radiative excitation at each level which results from the 
sum of radiation from the earth and lower atmosphere, the sun ana other 
levels in the atmosphere. 

The band functions are used to determine the effects o° 
absorption between a molecule at a given altitude and the various sources 
of radiation affecting a yiven band in the following way. The vertical 
optical paths to the lower boundary and to each level in the atmosphere are 
computed and, for daytime conditions, the optical path to the assumed top 
of the atmosphere in the direction of the sun using the secant of the solar 
zenith angle as a multiple. The effective radiation from the lower boundary 
is determined by multiplying the surface radiancy by the appropriate Ne func- 
tion. The effective radiation from the sun is determined by use of the 
proper "? function. The effective radiation from the other levels of the 
atmosphere may be obtained by integrating the product of emission at each 
level and the band function Neos 


tT 
= = 0 - 
Satm'*p) -f J(t, t) Noy (a, ‘) t) dt 
0 


. (43) 


b 
f J(t - ty) Noy (a, t - ty) dt 
p 
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Numerical integration is avoided for contributions from nearby 
atmospheric layers by assuming that the relative vibrational population between 
two points varies linearly with optical thickness. The integra’ relations be- 
tween the band functions are used in the form 


tT 


2 b 2 
f J(t) Noy (a, t) dt [ (A + Bt) Noy (a, t) dt 
ia 1 


= J(t,) yb (a, t,) - Ut.) nb (a, ta) 
V "12 1 2° "12 é 


J(t,) = (ty) 
+ at [us (a, tL) = hs (a, | (44) 


At large optical thicknesses, where the contribution to excitation is small, 
the radiative transfer functions used in the programs are not accurate enough 
for the subtractions in Equation (44) to yield meaningful results and Equa- 
tion (43) is evaluated by the trapezoidal rule. 

Vibrational populations are determined by balancing collisional 
and radiative excitation and de-excitation rates. Total band intensities 
at an observation point inside or outside the atmosphere are computed by de- 
termining the optical path from each level to the observation point HR 
integrating the product of local emission rate and the band function MP (a, Tt). 
This integration also uses a semi-analytical form for atmospheric layers close 
to the assumed detector location with the relation 


To 9 
[ Jt) MP (a, t) dt -[ (A + Bt) Mi) (a, t) dt 
1 TY) 


= Or) LP (ay xy) = ty) LP (a, 9) 


J(19) my J(r)) 
- Ty 


es (a, 1) - ys (a, 2)| (45) 
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where 


T 
LP (a, 1) “f ua, t) at (46) 
0 
Again, trapezoidal integration is used for layers far away from the detector. 

The spectral distribution of radiation is determined by 
finding a mean effective temperature and line-of-sight column density for 
the radiating molecules, computing intensity for each line in the band with 
the line function L(a. t), summing the line emission rates and scaling each 
line emission rate so that the total band intensity is obtained. 

The above analysis has been presented for the general Voigt 
profile. Actual use of the radiative transfer functions in the programs has 
been restricted to the Doppler profile. This restriction i. adequate for the 
water vapor and ozone bands which in the concentrations used in the program 
have a maximum vertical optical thickness at line center less than three. 

For limb viewing of carbon dioxide bands the Voigt profile should be used. 
Efforts to find a genera? way of setting up the radiative transfer functions 
to use a Voigt profile that varies with altitude have not been successful, 
Exploratory calculations show that the computed radiance values for strongest 
bands at a tangent height of 60 kilometers may be in error by 20 - 5U per- 
cent. In view of the uncertainties in chemical abundances, this error is 

not felt to be serious. 

The same radiative transfer functions have been used for 
the non-linear molecules as were used for the linear molecules. The match 
between radiative transfer functions has been made by determining a single 
"effective" rotational constant for the non-linear tri-atomic molecules. This 
was done by computing actual radiative transfer functions for a limited range 
of optical thicknesses adequate to cover that needed in the model. The com- 
puted functions were compared with those for the linear molecules and the 
effective rotational constant needed to produce an approximate match at optical 
thickness at line center of about two was used. 
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3.0 ATMOSPHERIC COMPOSITION 


There have been relatively few measurements of the variation with 
altitude of infrared emitting species at altitudes above 40 km, the limit 
of direct sampling by balloons. While ozone measurements by different 
worke’s agree reasonably well with each other and with theoretical computa- 
tions at altitudes up to 100 km, and nitric oxide profiles appear reasonably 
well known between about 80 and 150 km at middle latitudes, corresponcing 
data are Jess abundant for other species which are important. There is in 
some cases, indirect data available which when combined with computational 
models can yield limits to concentrations. Thus, Anderson (19712) has been 
able to put an upper limit on total hydroxy] radical abundance above 60 km 
from measurement of resonance fluorescence and with this result computes 
possible water vapor profiles between 50 and 100 km using a photochemical 
model that includes vertical transport by eddy diffusion. In a later paper, 
Anderson (1971 b) gives a measured altitude profile for hydroxy] with error 
limits + 120 percent. 

At present, then, estimates of chemical abundances of minor species 
recuired to predict infrared radiation require the use of photochemical 
computations. Even with the inclusion of best estimates of the effects of 
such transport processes as eddy mixing, these computations must be used 
with caution, even when they appear to agree with the limited concentration 
determinations that are available. Wofsy and McElroy (1973) give a discus- 
sion of methane concentrations in the stratosphere which illustrates the 
interplay of experiment and modelling studies, showing how estimates of 
eddy mixing coefficients may be obtained. 

While previous work (Degges, 1972, 1974) made use of a chemistry 
program which included one-dimensional diffusive transport, the current 
contract relies on published calculations to obtain estimates of the abun- 
dance of infrared emitting species in the atmosphere. The U.S. Standard 
Atmosphere, 1976 has been used as the source of abundances of the major 
atmospheric species No» 05 and 0. We discuss H,0, NO, 0, and CO, 
abundances below. 

Water vapor abundances are not well known, even in the stratosphere. 
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Balloon measurements are limited to altitudes below 40 km and are subject 

to considerable uncertainty as to the amount of water vapor carried to 

high altitude by the balloon. Direct measurements put lower and upper limits 
of the volume mixing ratio at two and ten parts per million in the upper 
stratosphere. A recent optical measurement from a balloon has been made 

by Gcldman, et al (1973) who observed thermal emission from the rotational 
spectrum near 25 micrometers. Near and below 30 km, the volume mixing ratio 
Js nearly three parts per million. In the region between 30 and 36 km, above 
the peak altitude reached by the balloon, the mixing ratio appears to be two 
or three times as large. Radford, et al (1977) have used grGund based micro- 
wave soundings to determine mesospheric water vapor concentrations. They 
obtain 10 ppm at 50 km with a maximum of 15 ppm at 60 km. Rogers, et al 
(1977a have analyzed radiance measurements obtained from the AFGL/DNA ICECAP 
Program to infer a water vapor mixing ratio of 3.5 ppmV between 50 and 70 km 
during late winter at Fairbanks, Alaska. They discuss other measurements and 
calculations. At present, we adopt 5 ppmV as a reasonable mixing ratio, as 
shown in Figure 4. 

A large number of measurements of nitric oxide concentrations are 
available. Representative mid-latitude concentrations are shown in Figure 5, 
which is based on the twilight measurements of Baker, et al (1977) and on the 
calculations of Oran, et al (1975). 

While ozone abundances below about 50 km must take into account 
lateral transport as well as vertical transport by eddy mixing and molecular 
diffusion, photochemical models are adequate above that altitude, and agree 
with observations obtained by measuring absorption of solar and stellar ultra- 
violet radiation and those deduced from measurements of the molecular oxygen 
1.27 wa infrared atmospheric band. The exact profile depends on assumed eddy 
diffusion coefficients, however, principally through atomic hydrogen and oxygen 
concentration. 0; abundances used here are shown in Fiqure 6. Below 60 kn, 
we use the ozone mode] in the U.S. Standard Atmosphere, 1976. These curves 
rely heavily on measurements by Llewellyn and Witt (1977) and Hays and Roble 
(1973). The calculations of Thomas and Bowman (1972) have been used as a 
guide to predict diurnal variations in abundances. Diffusive equilibrium is 
assumed above 100 km. 
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Figure 4: Water vapor concentration profile 
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Figure 5: Nitric oxide concentration values 
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Carbon cioxido concentration profile 
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Carbon dioxide has a lifetime of about ten days at the top of the 
atmosphere, where it is dissociated by solar radiation. It is well shicided 
by molecular oxygen betow about 120 km, however, and should remain well mixed 
below 100 km. The only good measurement above 100 kin remains that of Offer- 
man and von Zahn (1971), who find that between 120 and 140 km the mixing ratio 
is not more than 50 percent less than the tropospheric value. Above 140 kn, 
molecular diffusion is rapid enough compared to photodissociation to maintain 
diffusive equilibrium. Figure 7 gives our current estimate of CO, concentrations. 
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4.0 RESULTS FROM COMPUTATIONAL MODEL 
4.1 Introduction 


The revised high altitude spectral radiance programs can 
produce spectral radiance profiles for wavelengths between 2.5 and 25 
micrometers or between 400 and 4000 wave numbers. Results are presented 
here based on estimated abundance profiles of the four species CO, H0, 
0, and NO. The temperature and pressure profiles used were taken from 
U.S. Standard Atmosphere, 1976. Upwelling terrestrial radiances are assumed 
to be those corresponding to a black body source at the ground level 
temperature. 


4.2 Limb-Viewing Spectral Radiances 


Representative limb viewing radiance profiles are shown in 
Figures 8 through 19. These are based on the concentration profiles of 
Se.tion 3. The temperatures and number densities used are those of the 
U.o. Standard Atmosphere, 1976. 

Comparison of day and night profiles for a given tangent 
height shows that the low energy (small wave number and long wavelength) 
band radiances do not vary appreciably during the day. The high altitude 
high energy bands, however, do vary in intensity by as much as two orders 
of magnitude as the result of resonant and fluorescent emission following 
absorption of solar radiation. This is s*own clearly in the C0, 4,3 micro~- 
meter bands and in the water vapor complex at 2.7 micrometers. 


4,3 Individual Band Radiances 


For some p.“ ,oses it may be more convenient to have avail- 
able the total radiance of an individual band rather than the spectrum of 
tne region where the band is located. This may aid in making allowances for 
species abundances different from those assumed in a computation or in 
assigning an experimentally observed radiance value to one species rather 
than another. 
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Figure 8: Night limb viewing spectral radiance, 400 to 4000 cm”, 60 km 
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Figure 10: Night limb viewing spectral radiance, 400 to 4000 cm”, 80 kin 


tangent height. 
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Figure 11: Day limb viewing spectral radiance, 400 to 4000 om”! 
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Figure 12: Night limb viewing spectral radiance, 400 to 4000 cm™!, 100 kin 
tangent height. 
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Figure 13: Day limb viewing spectral radiance, 400 to 4000 on7!, 100 km 
tangent height. 
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Figure 14: Night limb viewing spectral radiance, 400 to 4000 on™! , 120 km 
tangent height. 


50 


: Pan es | 
SPECTRAL RADIANCE, WATTS cM”? ster™? (cm7?) 


tin a STR INN ee wae ee SU TBE ROT ater ARERR, NEE ee =~ ell ma tS Seip ah tee tain Ndi tne oa nena scent eae AA SBS RGR Fe Aine AE hE ace Anigae MMSE ae hg 


ae, anne am 


“12 


—_ 
So 


' 
— 
ww 


~ 
Qo 
‘ 
_— 
B - 


rr" 
So 
8 
os 
nn 


pwd 
oS 
t 
-_ 
an 


—_— 

27 
-_ 
~ 


10 
200 1200 2200 3200 4200 


WAVE NUMBER, CM! 


Figure 15: Day limb viewing spectral radiance, 400 to 4000 on™!, 120 km 
tangent height. 
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Figure 17: Day limb viewing spectral radiance, 400 to 4000 cm™!, 150 km 
tangent height. 
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Figure 18: Night limb viewing spectral radiance, 400 to 4000 cm), 200 km 
tangent height. 
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Figures 20 and 21 show the night and day total band radi- 
ances for the ten water vapor bands included in the spectra of Figures 8 
through 13. The 6.3 micrometer Vo band is the strongest, day and night. 
During the day, the V3 band at 2.66 micrometers strongly absorbs solar radi- 
ation and becomes almost as bright at a tangent height near 90 km. At night 
this band is nearly feur orders of magnitude fainter. 

Figures 22 and 23 show night and day limb radiance values 
for seven bands of the C0, 15 micrometer band compiex, the 4.3 micrometer 
Vy band and the two fluorescent bands originating from the V3 state, the 
10.42 micrometer 00011-10001 band and the 9.60 micrometer 00011. The 15 
micrometer band radiances vary little between day and night. During the 
day, however, the radiance of the 4.3 micromater band rises from its night 
value of between 1072 and 107° watts/cm@-ster to a value of between 10” 
and 2x 1077 watts/cm@-ster. This band is thick in the limb viewing 
geometry, and most of the increase is due to increased vibrational temper- 
atures at altitudes above 90 km. The 10.41 and 9.66 micrcemeter bands, which 
have much lower optical thickness, give a better indication of the alti- 
tudes at which the solar excitation is important. The 60 km value is about 
4x 1078 watts/cm@-ster day or night. Figures 24 through 27 show night and 
day band radiances for nitric oxide and ozone. 


4.4 Comparison With Observations 


A limited amount of data from the ICECAP Program is avail- 
able for comparison with the results of the radiance computer program. 
Figure 28 compares peak spectral radiance values for the 15 micrometer 
band reported by Stair et al (1974, 1975) and Rogers et al (1976, 1977b) 
with those computed by the present model and concentrations of Section 3 
(Model I) and the previous Model (Model II) using the U.S. Standard Atmo- 
spherc, 1962 temperature profile with 180 K mesopause temperature and uniform 
mixing of CO, at all altitudes. The differences between the models is due 
principally to differing temperature profiles and assumed CO, mixing ratios. 
A part of the difference between the data and the model can be resolved by 
using a temperature profile obtained at the observation site near the times 
of the flights. Complete resolution of the differences will require an 
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Figure 21: Day limb viewing total radiance values for water vapor bands. 


reat Pas 
i: One nee 


-_ 

‘x 

7 
y 

& 

Ee 

3 

3 14,98 wm 

g 

= 13.87 ym 


100 
TANGENT HEIGHT, KM 


Figure 22: Night limb viewing total raaiance values for carbon dioxide bands. 


59 


Spel 


ba Sea ae al Be Sle ek SR 


sits < 
ea re be Wteemamae ate soto 


* 


AO 


ETSI EN. 


LIM RADIANCE, WATTS #2 ster! 


16.74 15.45 


50 100 150 200 
TANGENT HEIGHT, KM 


Figure 23: Day limb viewing total radiance values for carbon dioxide bands. 
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Figure 25: Day limb viewing total radiance values for nitric oxide bands. 
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Figure 26: Night limb viewing total radiance value for ozone. 
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Figure 27: Day limb viewing total radiance value for ozone. 


improved model which takes into account the variation of line width and 
band shape with temperature in the region of the mesopause. 

Figure 29 compares the radiance model prediction for the 
ozone 9.6 micrometer band with data from the same sources as for CO,. The 
amount of ozone used in the calculation is relatively optically thin in 
the upward direction, and doubling the concentrations will approximately 
double the computed radiance. The model calculation agrees with observa- 
tion to within a factor of three for the 1973 measurement, and can be fit 
quite well simply by increasing concentrations above 80 km by a factor of 
three, with a smooth increase in concentration between 50 and 80 km. The 
shape of the observed radiance curve supports the qualitative features of 
the measurement by Hays and Roble (1973) and their interpretation (Roble 
and Hays, 1974) that the observed ultraviolet absorption is due to ozone 
and not a mesopheric aeroso] layer. 


4.5 Model Evaluation 


The most important general area of investigation in this 
research effort has been the improvement of models for the transport of 
infrared radiation. 

Figures 28 and 29 demonstrate that the radiance program 
results are in fair agreement with a small amount of observational data. 
The agreement with the carbon dioxide radiance model of James and Kumer 
(1973), is believed to be satisfactory. The model in its present state 
must be considered incomplete, however. For example, space and running 
time restrict the number of bands included, and further work is required 
to select a reasonable set of carbon dioxide isotopic bands for inclusion. 

The program results are highly dependent on the atmospheric 
temperature and concentrations used, as is evident from Figure 28. Uncer- 
tainties in atmospheric composition, which commonly range from factors of 
two to orders of magnitude, are reflected in the computational results, 
and produce the greatest uncertainty at high altitudes. Uncertainties in 
the “ates of collisional excitation mechanisms also affect the calculations. 

The program for computing band radiances is now slightly 
faster than the 1974 mode. The most time consuming process is now the 
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integration of band radiances. The spectrum generated by the program in 
its current form does not properly represent line intensities for bands 
where a warm emitting region is overlaid by 2 cooler abserbing region. 
Future work will include use of the FASCODE algorithm, described in Section 
6, to improve the calculation of both the radiative transfer near the meso- 
puase and the spectral calculations. As in the past, results of future 
radiance measurements will guide in program modification. 


ok te 


5.0 VIBRATIONAL TEMPERATURES 


The recognition that vibrational temperatures at high altitude are 
not the same as the local translational temperatures is one of the most im- 
portant reasons for the development of radiance models such as this one and 
that of James and Kumer (1973). The present program has been adapted to 
provide vibrational populations and temperatures at altitudes between ground 
level and 700 km. No one set of radiative transfer functions is adequate 
to compute outgoing radiances for the full range of altitudes, but even the 
Doppler profile can be used to estimate departures of vibrational tempera- 
tures from kinetic temperatures. 

Figures 30 through 33 present vibrational temperatures for one or 
two of the lowest excited vibrational states of carbon dioxide, water vapor, 
ozone and nitric oxide. 

Figure 30 shows day and night vibrational temperatures for the 
01101 Vo mode excited level and for the 00011 V3 mode excited levels. The 
Vv, mode remains close to thermal equilibrium up to about 75 km altitude. 

The V3 mode temperature remains close to that of molecular nitrogen at alti- 
tudes up to about 70 km. Kumer and James (1974) have shown that the nitro- 
gen and the V3 mode vibrational temperatures depart from kinetic temperatures 
above 50 km, while remaining closely coupled to each other. The 10001 aT 
mode level (not shown) is not in thermai equilibrium above 30 km, and higher 
vibrational levels depart from thermal equilibrium faster than the lower 
levels. 

Figure 31 shows temperature for the water vapor Vo mode with radi- 
ation near 6.3 micrometers. Thermal equilibrium is maintained to near 70 
km. Resonant absorption of sunlight raises the daytime vibration temperature 
a few degrees above the nighttime temperature above 70 km. 

Figure 32 shows temperatures for the ozone V3 mode, which remains 
in thermal equilibrium up to 60 km. The peak in the daytime temperature 
comes from the chemiluminescent formation of ozone by the three body reaction 
of atomic and molecular oxygen. The rate constant used is that of Garvin 
and Hampson (1974). 
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Figure 31: Water vapor Vo mde vibrational temperatures 
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Figure 33: Nitric oxide vibrational temperatures 
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Figure 33 shows nitric oxide vibrational temperatures. Vibrational 
levels are not in thermal equilibrium above about 10 km. The temperatures 
between 10 and 30 km may not be correct. The concentrations used in the 
calculation are only estimates and no attempt was made to differentiate 
between day and night concentrations at low levels. The vibrational excita- 
tion rates used in the calculations are taken fron Taylor's (1974, 1975} 
compilations, except for the use of Gl¥nzer and Troe's (1975) and Quack and 
Troe's (1975) rates for excitation of nitric oxide by atomic oxygen. The 
results using these rates appear to confirm the suggestion made by Degges 
(1971) that atomic oxygen should be an important collision partner above 
the mesopause. 

A related computer program (Hurd, et al, 1977) has been developed 
to follow the time history of vibrational populations of the C0, V3 mode 
(the 0011 level) and molecular nitrogen and oxygen vibrational temperatures. 
This gives more accurate populations than the BGND Program can because BGND 
must assume steady state conditions and does not allow for the possibility 
of nitrogen vibrational temperature variation. In addition to its time 
history capability, use of matrix solution of the radiative transfer problem 
has been investigated using NORAD, with the conclusion that execution of the 
BGND Program will be more efficient using nitrogen vibrational populations 
computed by NORAD. 
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6.0 FASCODE — A LINE-BY-LINE ATMOSPHERIC TRANSMISSION AND 
EMISSION COMPUTER PROGRAM 


A line-by-line computer program, FASCODE, for the calculation of 
atmospheric transmission and emission has been developed. This has in- 
volved the extension of the HIRACC algorithm developed at AFGL for the 
rapid line-by-line calculation of atmospheric infrared spectral absorption 
coefficients at a given temperature, pressure and composition assuming a 
Lorentz shape for the spectral line profile. A complete description and docu- 
mentation will be found in an accompanying scientific report, Smith et al (1978). 

The Lorentz shape is not appropriate at higher altitudes where 
Doppler broadening becomes dominant. To enable the use of the code at all 
altitudes a Voigt profile version of the HIRACC algorithm wes developed. 

This program yields a line shape smoothly varying from nearly Lorentz at 
lower altitudes to purely Doppler at altitudes 2 50 kms. In addition, the 
domain of the line shape function was extended from 48 to 64 half-widths. 
With these new features the extended HIRACC algorithm required only a 
moderate increase in running time over the Lorentz profile version. 

Using a layered atmospheric model, the spectral absorption prop- 
erties of each layer are different. To compute spectral absorption coef- 
ficients over a given path through these layers requires a careful treat- 
ment since the line widths and/or shape for each spectral line will, in 
general, change with altitude. The HIRACC algorithm applied at two dif- 
ferent altitudes will produce results at different spectral sampling intervals 
requiring an appropriate method of merging these results in order to produce 
an absorption coefficient for the given path. Since one needs a coarser 
sampling interval at the lowest altitudes of the sight path chosen, the 
merger algorithm is started at this point. The HIRACC algorithm is exercised 
and the results are written to mass storage. Using criteria established for 
sampling intervals which can be handled by the merging algorithms, the 
resolution for the next level is chosen and the HIRACC algorithm is 
exercised. The results at this level are then merged with the previous 
level (at a coarser sampling interval) by interpolating the previous results 
into the new. This procedure is followed until ali levels have been completed. 
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The resulting spectral absorption signature is then available for use for 
the particular path calculated. For convenience a plotting program has 
been included within the FASCODE structure. We note here that the merged 
absorption coefficient obtained is valid for a path proceeding (“Tooking") 
in either direction. 

The emission along a given path may also be calculated if one 
assumes thermodynamic equilibrium along a path. Using this assumption an 
option has been included in FASCODE for the calculation of emission as well 
as transmission. Emission calculations require a new feature however, since 
emission along a path looking up is not in general equal to the emission 
along the same path looking down. In order to reduce calculation time, the 
emission calculation proceeds from the lowest level (coarsest spectral 
interval) upward to the highest level (finest spectral interval) regardless 
of the direction of sight along the path, but the algorithm for the emis- 
sion calculation is different for the upward and downward looking cases. 
Merging of the results for the layers is done in the same manner as for 
the absorption, except that for this case the transmission is calculated 
instead of the absorption and both the transmission and the emission are 
produced. 

The results of the programs have been tested for internal self- 
consistency and compared with results from other line-by-line calculations, 
in particular the work of Kyle and Goldman (1975). Satisfactory agreement 
has been obtained in all cases tested. The interested reader is referred 
to the detailed scientific report describing the method, Smith et al (1978). 
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7.0 HIGH ALTITUDE RADIANCE PROGRAMS 


This section lists and describes the main programs revised ur 
written for this contract, BCKGND and SPCTRA. 

BCKGND computes vibrational populations and total band radiances 
for an earth's limb viewing geometry at 2 kn intervals between 0 and 150 
and at 10 km intervals between 150 and 700 km. 

SPCTRA computes spectral radiances corresponding to the total 
band radiances output by BCKGND. Control information is read from TAPE 5 
= INPUT. Permanent data is read from TAPE 1 and the results of the BCKGND 
computation are read from TAPE 3. Qutput data is sent to TAPE 6 = OUTPUT 
for printing and to TAPE 4 for use by a plotting program. Optional output 
of unf?)tered spectral radiances may be placed on TAPE 7 and saved for 
later use. 

Listings of two auxiliary programs, INTERP and NORAD are included 
for information, but are not documented further. 

INTERP produces molecular concentrations for No» 0, and 0 consis- 
tent with the U.S. Standard Atmosphere, 1976. The ozone concentrations of 
Figure 6 are included in DATA statements for arrays O3DAY and O3NITE in Sub- 
routine 03. The original version of INTERP was used by Corbin, et al (1969) 
to generate mode] atmospheres for low altitude limb radiance calculations. 
Part of the output from INTERP is used as input to Program BCKGND. The sub- 
routine with the same name is no longer needed in BCKGND. 

NORAD is used to compute time dependent vibrational populations of 
No and the C0, V3 mode. Its use in auroral calculations is illustrated in 
Hurd, et al (1977). 


17 


PROGRAM BCKGND (BGND) 


Purpose: 


Program Flow: 


The main Program BCKGND (short for BACKGROUND) serves pri- 
marily to control the operation of subroutines which com- 
pute vibrational populations and radiative rates. BCKGND 
reads and checks the input data, and transfers contro] to 
Subroutine BANRAD which computes vibrational populations, 
local optical thicknesses and volume radiances for all 
included bands of a single molecular species. When con- 
trol is returned, BCKGND computes and outputs integrated 
band radiances for ‘imb viewing (exo-atmospheric) and 
seven endo-atmospheric viewing angles at seiscted tangent 
heights or endo-atmospheric altitudes. 


The initial step in operation of BCKGND is the reading of 
all input data. The number IGAS of molecular species 
including total number density is read in. Then in the 
DO range ending with Statement 170, the following opera- 
tions are performed. The title card, containing alpha- 
numeric identification, an index check number and maxi- 
mum altitude considered for the species is read in and 
its check number checked against the index of the DO 
range. If the check succeeds, altitude and number den- 
sities are read into arrays ALT and STORE until an 
altitude corresponding to the maximum altitude is found. 
hen this altitude is found, the logarithm of input number 
densities is formed in the DO range ending with Statement 
130 and contro] is passed to Subroutine INTERP which 
interpolates to form values at 1 km intervals up to 160 
km and at 5 km intervals at higher altitudes. The inter- 
polated values are converted to number densities in the 
DO range ending with Statement 160. In the DO range 
ending with Statement 175, total atmosphere number densi- 
ties are stored in array TOTDEN. Temperature data are 
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then read in similarly in the DO range ending with State- 
ment 245 and interpolated values found by Subroutine INTERP. 
The complete set of number densities and temperatures is 
then printed out. Next, vibrational level and band con- 
stants are read in for each species in the DO range ending 
with Statement 331. The first card for each species con- 
tains an alphanumeric identification, number of vibrational 
levels NLEV, number of bands NBAN, molecular mass AMAS, 

and two vibrational energy transfer coefficients DEX and 
DEXB. The alphabetic code is checked and if it is a valid 
code stored in array DATATY, these constants are transferred 
to the appropriate arrays. In the DO range ending with 
Statement 313, constants for each vibrational level are 
read in, and card sequence numbers checked. In the DO 
range ending with Statement 316, constants for each band 
are read in and card sequence numbers checked. For C0, 

and H0, tabular values of temperatures and associated 
vibrational energy exchange coefficients are read before 
the level and molecular band data. 


The computation uf molecular band radiances then proceeds 
in the following manner. The outer DO range index NLTE 
and ending with Statement 900 controls the factors con- 
sidered in a complete set of calculations for any of five 
possible sets. The corresponding DO statement may be 
varied depending on what factors are desired. See the 
description of Subroutine BANRAD for the effect of dif- 
ferent values of the index NLTE. Each cycle through the 
nested DO range with Index I and ending with Statement 
850 results in the computation of band radiance values 
for a single molecule. Action within this DO range is 

as follows. 
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Subroutine BANRAD is called and returns volume emission 
radiance values in array RAD and local optical thickness 
values in array TAU. Then for each band, line-of-sight 
band radiances are computed in the DO range ending with 
Statement 840. Within this DO range, the value of NLTE 
is examined and appropriate headings are printed. Then 
trapezoidal rule integration is performed to determine 
the radiance arising along each of the eight lines-of- 
sight, including the modification necessary for the 
effects of optical thickness along the line-of-sight. 

The DO range ending with Statement 830 controls the inte- 
gration for each tangent height and endo-atmospheric 
viewing level. In the current version, radiance values 
are computed at 5 km intervals of tangent height up to 
150 km and at 25 km intervals above this level. The 
integration of radiance values is performed for each 
endo-atmospheric viewing level in the DO range ending 
with Statement 810. The integration begins at the view- 
ing altitude. At each point on the line-of-signt, the 
contribution of the local volume emission rate to that 
received at the viewing point is computed by Subroutine 
SMI, which takes into account the effects of optical 
thickness between the emitting point and the viewing 
point. The equivalent integration for the limb viewing 
case is performed in the DO ranges ending with Statements 
815 and 818. Upon completion of the integration for each 
tangent height and endo-atmospheric viewing level, band 
radiances for that level! are punched and printed. This 
punched output, together with that produced by Subroutine 
BANRAD, is used as input data to the special radiance 
Program SPECTRA. 
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Tables 7.1 through 7.4 list the molecular bands and 
related physical data now included in the programs. 
Band strengths for water vapor, carbon dioxide and 
ozone are those used in the AFGL absorption line 

atlas described by McClatchey, et al (1973). Band 
strengths for nitric caide are approximately those 

of Abels and Shaw (1966) and Schurin and Ellis (1966). 
Infrared solar fluxes at band centers are interpolated 
from the tabulation of Johnson (1965). The tables are 
followed by descriptions of the subroutines, the Pro- 
gram Listing and Sample Input Data. 
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TABLE 7.1: WATER VAPOR BAND DATA 


Vibrational 


Band Strength Solar Flux 

jue oa iavelenath _ ae K * Band oer 1 
State State (mol ~ cm" cm) (photons cm“ pm ~) 

010 000 6.27 1.06(-17)* 5.60(15) 

020 000 3.17 6.58(-20) 2,60(16) 

100 000 2.73 3.62(-19) 5.10(16) 

001 000 2.66 7.99(-18) 5.40(16) 

020 010 6.42 8.61(-21) 5.70(15) 

100 010 4,85 8.90(-23) 1.20(16) 

001 101 4.63 4,90(-22) 1.50(16) 

030 010 3.26 7.99(-23) 3.30(16) 

110 010 2.75 1.50(-22) §.00(16) 

011 010 2.68 3.30(-21) 5.40(16) 

030 020 6.60 7.00(-24) 4,90(15) 

030 000 2.14 2.00(-22) 9.30(16) 

110 000 1.91 1.83(-20) 1.20(17) 

011 000 1.88 9.16(-19) 1.25(17) 


*Number in parentheses is power of ten. 
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TABLE 7.2: CARBON DIOXIDE BAND DATA 


Vibrational 


Tete (Poe hss aN eG, So BT ay age te et Sek See 


: Band Strength Solar Flux 
umper” Cover Wavetempey, «at 295.K Band Center 
State State (mol ~ cm” cm”) (photors cm ~ pm) 
01101 00001 14,98 8.26(-18) 4.60(14) 
10002 01101 16.18 1.44(-19) 3.70(14) 
02201 01101 14.98 6.49(-19) 4,60(14) 
11102 10002 15.45 2.22(-20) 4.20(14) 
11102 02201 16.74 5.21(-21) 3.35(14) 
03301 02201 14.97 3.82(-20) 4.60(14) 
10001 01101 13.87 1.85(-19) §.80(14) 
00011 10001 10.41 4.91(~22) 1.31(15) 
00011 00001 4.26 9.60(-17) 1.60(16) 
00011 10002 9.40 6.32(-22) 1.78(15) 
10012 00001 2.77 1.04(-18) 4,.95(16) 
10011 00001 2.69 1.50(-18) 5.35(16) 
10012 10002 4.30 1,93(-19) 1.60(16) 
10011 10002 4.12 1.06(-22) 1.75(16) 
10012 100u1 4.50 1.27(-22) 1.39(16) 
10011 10001 4.30 1.18(-19) 1.60(16) 
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Vibrational 
Transition 
Upper Lower 
Stat State 
1 0 
2 0 
2 1 
Viorational 
Transition 
Upper Lower 
state State 
001 000 


TABLE 7.3: 


Band Center me tet dah 

Wavelength PS ee ees | 

(ya) (mo on® om") 
5.32 5.0 x 1078 
2.68 8.8 » 10°20 
5.41 1.27 x 1072! 


TABLE 7.4: OZONE BAND DATA 


Band Strength 
Band Center 
Wavelength ene “1 
—m) ss (mol cn” cn") 
9.60 1.29(-17) 
354 


NITRIC OXIDE BAND DATA 


Solar Flux 
Band Center 


(photons on”? ym”) 
8.7 x 1028 
5.4.x 1086 
8.4 x 108 


Solar Flux 
Bend Center 


(photons om”? yn” ') 


a eke we 


FLOWCHART: 8CKGND 


SFT CONSTANTS FOR INPUT AND CUTPUT CHANNELS 


READ TEMPERATURE AND MOLECULAR CONCENTRATION DATA | 
READ CONCENTRATION AND TEMPERATURE DESCRIPTORS 


LIST CONCENTRATIONS AND TEMPERATURES | 


BEGIN DO 331 LOOP TO READ PHYSICAL DATA FOR ALL MOLECULES 


we 


| DO 213 LOOP TO READ CONSTANTS FOR EACH VIBRATIONAL LEVEL 
| DO 316 LOOP TO READ CONSTANTS FOR EACH INFRARED BAND 


END 00 331 LOOP 


SET INITIAL VALUE OF NLTE FOR TYPES OF CALCULATIONS 


SET INITIAL VALUE OF I FOR DESIRED SPECIES 


0 Seeppmenennmse 


[WRITE L LEVEL AND BAND DATA FOR INPUT TO SPECTRA | 


= 


uF 


ENTER SUBROUTINE BANRAD TO COMPUTE VIBRATIONAL POPULATIONS 


SET J=1 TO BEGIN CALCULATIONS FOR INDIVIDUAL BANDS 


PRINT PAGE HEADINGS 


DO LOOP 830 TO INTEGRATE RADIANCES FOR EACH 
TANGEN HEIGHT AND 7 ENDOATMOSPHERIC ANGLES 


DO LOOP 818 TO INTEGRATE RADIANCES FOR LIMB VIEWING GEOMETRY | 


PRINT BAND RADIANCES 


| WRITE BAND RADIANCES FOR INPUT TO PROGRAM SPCTRA 


Subroutine BANRAD 


Purpose: 


Method: 


BANRAD computes the population of vibrational levels of 
infrared active species, including the effects of collis- 
jonal excitation; absorption of radiation from the lower 
atmosphere, sun and molecules radiating at high altitude. 
The values returned to the main program are volume radi- 
ance and unit optical thickness at each altitude for each 
band included in the computation. 


Subroutine BANRAD computes the population of vibrational 
levels of infrared active species, including the effects 
of collisional excitation and absorption radiation from 
the lower atmosphere, sun and molecules radiating at high 
altitude. The values returned to the main program are 
volume radiance and unit optical thickness at each alti- 
tude for each band included in the computation. 


The collisional processes which are assumed to effect the 
vibrational state of a molecule are (1) translational- 
vibrational (T-V) interactions in which translation energy 
is transformed into vibrational energy or vice-versa dur- 
ing @ collision, (2) intermolecular vibrational exchange 
(V-V) in which a molecule of one species exchanges a quan- 
tum of vibrational energy with a molecule of another species, 
and (3) intramolecular vibrational exchange of polyatomic 
molecules in which a collision transfers the energy in one 
mode of vibration into another mode within the same 
molecule. 


When a molecule in the ground state collides with another 
molecule there is a finite probability Pry that a portion 
of the relative translational energy will be transformed 

te vibrational energy if the translational energy is equal 
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to or exceeds the vibrational energy. Similarly, if the 
molecule is in an excited vibrational state, there is a 
finite probability Pur that the vibrational energy will be 
converted to relative translational energy during the col- 
lision. For a Boltzmann distribution of translational 
energies at a given temperature 


¢ 9y he 
Pry * Pur" g, exe (: str 


where 9, and 9, are the statistical weights of the upper 
and lower vibrational levels, h is Planck's constant 
(6.6256 x 10°27 erg-sec), c is the velocity of light 
(2.9979 x 1910 cm/sec), k is the Boltzmann constant 
(4.28054 x 19°46 erg/K), T is absolute temperature (K) 
and A, is the wavelength (cm) corresponding to the 

energy difference of the two vibrational levels. Similar 
expressions hold for intermolecular and intramolecular 
V-V transfer. 


The rate at which a vibrational level of a single molecule 
is excited (Kry) or de-excited (Kyy) is given by the pro- 
duct of the collision frequency Z and the transition prob- 
ability. The subroutine uses temperature dependent rate 
coefficients Krys kyr etc., such that if [M] is the con- 
centration of the collision partner, Kyy = ky CM], etc. 


The radiative processes which are assumed to affect the 
vibrational state of a molecule are (1) adsorption of solar 
radiation, (2) absorption of radiation from the lower atmo- 
sphere, (3) absorption cf radiation emitted by molecules 

in the upper atmosphere and (4) spontaneous emission at a 

rate A (sec ”), the Einstein A coefficient. Because ab- 
sorbing molecules are present between each source of radiation 
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and the molecule for which the vibrational state is to be 
computed, radiative transfer functions must be used to take 
into account the effects of absorption. The mode] assumes 
that individual lines have a Doppler (thermal velocity 
broadened) spectral contour. 


At a given altitude, there will be Ny molecules/cm* in the 
ground vibrational state of a molecule and Nu molecules/cm> 
in an upper state. Molecular nitrogen and oxygen are assumed 
to be the only collision partners effective in transfer of 
vibrational quanta, and their concentrations are used to de- 
termine total vibrational excitation rates Kry> Kyy> etc. 

The treatment below refers to a two-state molecule. 


In unit volume, molecules in the ground state are excited 
to the upper state by: 


(1) Collisional T-V excitation at a rate Nikry 


(2) Collisional V-V excitation at a rate Ny Kay 
(the subscript MV is used for excitation, 
VM for de-excitation, by transfer of vibra- 
tional quanta). 
(3) Absorption of solar flux at a rote NyI,G'M? 
(k,) where G' is molecular cross section for 
absorption of radiative flux. 


(4) Absorption of lower atmospheric radiance at 
b 
a rate NaN, eG 'Ny o(ky) 
(5) Absorption of radiance emitted at other levels 
of the atmosphere at a rate Nf Ny 16'No, (kp) dz 


Molecules in the upper state are de-excited by: 
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(1) Collisional V-T de-excitation at a rate 
Nukvr 

(2) Collisional V-V de-excitation at a rate 
Nuke 


(3) Spontaneous radiative de-excitation at a 
rate Nua 


Excitation rates are assumed to balance de-excitation rates 
and the equation 


ees ae + 1.6'M(k en GN. (k.) 
7 a (ery Kay + TyG'My Che) + Ny gSONy tke 


b 
+ fy BND (kg )82 - Ny (Kyy + Ky ¢ An) =o 


1s solved by iteration to determine Ny and Ny at all 
altitudes. 


thermal equilibrium populations of vibrational levels. This 
{s performed in the DO range ending with Statement 130. 

Each passage through tiris D0 range computes populations for 
one altitude. Within this DO range are two cthers. The 

DO range ending with Statement 110 computes the vibrational 
population of each level relative to the ground vibrationai 
population of each level relative to the ground vibrational 
level, using the energies (wave number) stored in WAVE and 
the statistical weights stored temporarily in real array 
VIBPOP. The sum of relative populations is stored in SUM. 
The final ground state population is then computed by divid- 
ing total concentration by SUM, and stored in real array 
GNOSTE. The DO range ending with Statement 120 then computes 


| Program Flow: The first step in the computation is determination of the 
| 
j 
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the vibrational pepulation of each level by multiplying total 
concentration by relative concentration divided by SUM. 


The 90 range ending with Statement 135 selects the vibra- 
tional level which is coupled with the first excited nitro- 
gen vibrational level. 


The 00 range ending with Statement 160 computes the rate 
at which a single molecule in a given level is excited 
or de-excited by collisions at each altitude. 


The outer DO range ending with Statement 180 computes wave- 
lengths (BNDLAM), Einstein A coefficients (EA), rates for 
absorption of radiation from the sun (GDOWN) and from the 
earth and lower atmosphere (GUP) and the G factor used in 
computing radiative transfer within the upper atmosphere 
GFAC). The inner DO range ending with Statement 180 com- 
putes the unit optical thicknesses for each band at each 
aii:tude using Function TAUMAX, and stores the result in 
array BNDTAU. In addition, it computes the volume radi- 
ances and stores them in array BNDRAD. 


At this point, if NLTE is 1, computation for thermal equi- 
librium is complete and control] is transferred to Statement 
900. If NLTE is not 1, the computation for non-equilibrium 
is performed iteratively in the 00 range ending with State- 
ment 400. The DO range ending with Statement 190 computes 
various integrated optical thicknesses for each band. The 
DO range ending Statement 182 computes the total thickness 
between each two altitudes in ther vertical direction and 
stores the values in array TTH. 


The DO range ending with Statement 184 (1) computes the 
optical thickness between the base altitude and each other 
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altitude, storing the result in array TTUP and (2) computes 
the total optical thickness between the top altitude and 
each other altitude, storing the results in array TTDOWN. 

If the computation is not for daytime conditions, the DO 
range ending with Statement 186 sets the values in TTDOWN 
equal to 0.0, for efficiency in the later computations. If 
the computation is for the effects of collisional excitation 
only, the DO range «nding with Statement 188 sets the values 
In TTUP equal to 0.0. 


The computation of vibrational populations now begins. 
Control is transferred to 

Statement 210 if I is 1, CHy 

Statenent 220 if I is 2, C0, 

Statement 230 if I is 3, H0 

Statement 240 if I is 4, NO 

Statement 250 if I is 5, No0 

Statement 230 if I is 6, pre-dawn Ozone 

Statement 230 if I is 7, noon Ozone 


Following each of these statements is a DO range in which 
vibrational populations are computed. A discussion of two 
statements in the computation for C0, will illustrate the 
computational details. First, look at the statement in the 
DO range immediately following Statement 221 that begins 
VIBPOP(NZ,8) = which computes the population of the first 
level of the V3 vibrational made. This vibrational level 

is coupled to the nitrogen first vibrational level by col- 
Tision, to the upward flux from the earth and lower atmo- 
sphere, to solar flux, and to lower vibrational levels by 
radiation. Excitation by absorption of radiation from within 
the upper atmosphere is neglected for this level. A ground 
state molecule is excited by collision with nitrogen at the 
rate VVNSIG(NZ) (collision exchange), by absorption of solar 
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flux at a rate GDOWN(9)*SM1(KB(9), TDOWN2(NZ)) (equal: 
1,6°E(k,)) and by absorption of flux from the earth and 
lower atmosphere at < rate GUP(9)*SLM2(KB(9), TTUP2(NZ)) 
equals Ny gS'ND(Ky))- An excited molecule is de-excited 

by radiation at the rate EA(9) (other Einstein A coefficients 
are negligible), by vibrational exchange with nitrogen at a 
rate VVRSIG(NZ), and by transfer of vibrational energy to 
other modes within the CO, molecule at a rate VVDSIG(NZ). 

The vibrational population is then the ground state popu- 
lation times the ratio of excitation and de-excitation rates. 
Second, look at the statement beginning POP2 = ..., which 
computes the vibrational population of Level 2, the first 
excited level of the Vo vibrational mode. This level is 
excited from the ground state by collisional excitation at 

a rate TVS = TVSIG, by absorption of radiation from the lower 
atmosphere and earth at a rate GUP(1)*SLM1(KB(1), TTUP(NZ1)) 
(equals Ny gS'ND(Ke)), by absorption of radiation from other 
points in the upper atmosphere in the 15.0 um fundamental 
band at a rate 2nGFAC(1)*(XFER(VIBPOP(1,2), TTH(1,1), NZ, 

1, 1, KB(1)) + XFER(VIBPOP(NZ,2), TTH(NZ,1), NALT, NZ, 2, 
KB(1)))*EA(1) (equals J ,p8"N2(kp)¢2) and by absorption 


of solar radiation at the rate GDOWN(1)*SMI(KB(1), TTDOWN 
(NZ,1)) (equals 16°, (k,)). It is excited from Levels 3, 

4 and 7 by collisional excitation VTS and by radiavion EA(2), 
EA(3) and EA(7). De-excitation occurs by radiation EA(1), 

and by collisional de-excitation to the ground state at a 

rate VTS = VTSIG(NZ), and by collisional excitation to higher 
states at a rate TVS. The approximate vibrational population 
is then the ratio of the sums of population rates times higher 
state populations divided by the sum of de-excitation rates 
for single molecule. 


After each set of vibrational populations is completed, new 
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values for volume bard radiances are computed in the nested 
DO ranges ending with Statement 360. After three iterations 
te compute vibrational populations, contro] is returned to 
the main program. 


inputs: NLEVEL, NBAND, NALT, CDEXT, CDEXB, RM, LSC, LCC, I, LBC, 
LBU, BV, STR, WAVE, CONC, SFLUX, TFLUX, KCHEM, NR, NP, NW, 
NLTE, TEMP, BNDLAM, DENTOT, NBAND1. 

Outputs: VIBPOP, BNOTAU, BNDRAD 

MNEMONIC VARIABLES 

AA intermediate factor in computing vibrational populations 

AAA quantity used in NO chemiluminescence calculation in original 
Subroutine BANRAD application; it has no effect at present 

BNOL wavelength of current interest (ym) 

BNOLAM wavelength array (ym) 

BNORAD volume radiance 

BNOTAU unit optical thickness 

B rotational constant (cn?) 

CDEXB rate coefficient at 300 K for vibrational de-excitation 
of lowest NO vibrational level (om? sec™!) 

CDEXT rate coefficient at 300 K for vibrational exchange of NO 
with No or 05 for the most significant vibrational level 
(em? sec”!) 

CHEML currently inactive 

CONC array of NO concentrations (em™?) 

C0, real array used to store reduced temperature value and rate 
coefficients to generate vibrational energy exchange rates 
for C0,: currently inactive 

DENTOT array of total number densities at altitude mesh (cm79) 

EA array of Einstein A coefficients 
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FAC 
FACTOR 
GDOWN 


GFAC 


GNOSTE 


GUP 


NB 
NBAND 
NBANC 1 
NBC 


used as temporary storage 

temporary storage in computing vibrational exchange 
real array in which are stored rates of absorption 

of solar radiation, 16" for each band 

real array in which are stored values of S' for each 
band 

real array in which are stored ground state vibrational 
populations, equivalenced to VIBPOP(1,1) 

real array in which are stored rates of absorption from 
the lower atmosphere, 2nN) 4G" for each Sand 

index of species of interest (equals 4 for NO) 

integer used as index in iterative computation of 
vibrational level population 

index for altitude level 

index for final vibrational level in band transition 
index for initial vibrational ‘evel in a band transition 
integer used to compute JI, JF 

index for vibrational level 

integer array used to store band types 

index of band type 

chemiluminescence index; currently has no effect 
integer array storing level types 

current level type 

integer array containing information on which levels pro- 
duce a band 

integer array containing band type information 

integer array containing vibrational level coupling go 
No or 05 

integer array containing level type information 

number of altitudes for which results are to be returned 
index for band number 

number of bands for this species 

equal to NBAND 

temporary storage index 
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NBX 
NF 

NI 

NL 
NLEVEL 
NLTE 


NMIN 
NP 

NR 

NW 

NY, NZ 


POP2, POP3, 
POP4, POPS, 
POPS, POP7, 
POPS 


POWER 


temporary storage index 

temporary storage index 

temporary storage index 

level number index 

number of vibrational levels for this species 

integer variable used to indicate conditions of computa- 
tion, 1 for local thermodynamic equilibrium, 2 for 
collisional excitation only, 3 for night conditions, 4 
for noon conditions, 5 to include high vibrational temper- 
ature for nitrogen 

integer variable storing (NALT - 1) 

output device index 

input device index 

output device index 

index for altitude level 

real variables used temperature storage for vibrational 
populations 


real array in which are stored conversion factors between 
emission rates in photons om”? sec”! and watts cm”? ster! 
for each band 

molecular mass for this molecule (g/mole) 

array of solar flux values for each band 

temporary storage of band strength (em™! (atm-cm) =) 

array of band strength for each band (em’ }(atm-cm)~!) 
used to store partial sums 

real array equivalenced to TTDOWN(1,9) 

array containing temperature profile (Kk) 

array containing fluxes for lower altitudes for each band 
temporary storage for temperature (K) 

nitrogen vibrational tempzrature (K) 

array storing values of total optical thickness from top 
of atmosphere for each band 
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TTH 
TTUP 
TiUP2 
TVS 
TVSIG 
VIEBPOP 
VTS 
VTSIG 
WVDSIG 


WNSIG 


VVRSIG 
WUSIG 


VVWAVE 
WATVA 
WAVE 
XYM 


XYMM 
XYN 


Non-Librar 
Externals: 


Called by: 


array storing values of total optical thickness between 
adjacent altitudes for each band 

array storing values of tctal optical thickness from 
bottom of upper atmosphere for each band 

array equivalenced to TTU: (9,1) 

used to store local value of TVSIG 

array storing values of vibrational populations of each 
state for each altitude 

array storing values of vibrational populations of each 
state for each altitude 

local value of VTSIG 

array of vibrational de-excitation rates at each altitude 
array of estimated values of transfer of vibrational 
energy (downward) 

array of rates of transfer of vibrational energy from 
No to this molecule 

array of rates of transfer of vibrational energy to No 
array of estimated values of rate of transfer of vibra- 
tional exchange (upward) 

difference between No vibrational energy and a molecular 
vibrational energy 

array storing reduced temperature values and rate coef- 
ficients to generate vibrational energy exchange rates 
for C0, 

array containing energy level of each band (cm™?) 
temporary storage variables 

temporary sturage variable 

temporary storage variable 


EVAL, SIMP, SLM1, SM1, TAUMAX, XFER 


BCKGND 
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START 


DO 130 LOOP TO INITIALIZE POPULATION 
ARRAY VIBPOP ASSUMING EQUILIBRIUI 


DO 135 LOOP TO INITIALIZE VIBRATIONAL TRANSFER VARIABLES | 


FLOWCHART: BANRAD 


DO 160 LOOP TO COMPUTE VIBRATIONAL TRANSFER anys | 
RSIS eee irra at rane 
170 


DO 180 LOOP OVER NUMBER OF BANDS TO COMPUTE 
RADIATION FACTORS, OPTICAL THICKNESSES AND BAND RADIANCES 


ITER=ITER+1 


| GO 190 LOOP TO COMPUTE ARRAYS OF TOTAL OPTICAL THICKNESSES 


 eeenanennedinahadaemamee’ 


COMPUTED GO TO BRANCH TO THE APPROPRIATE CALCULATION 
OF THE ARRAY OF VIBRATIGNAL POPULATIONS INCLUDING 
EXCITATION AND DE-EXCITATION MECHANISMS FOR EACH SPECIES 


DO 360 LOOP OVER NUMBER OF BANDS 10 COMPUTE 
OPTICAL THICKNESSES AND RADIANCES 


tones 


Function XFER 


Method: 


Program Flow: 


Function XFER is used to compute the excitation of a vi- 
brational level caused py "olecules radiating in a given 
band. 


Function XFER computes the contribution to the source func- 
tion from the radiation received in 4 given vibration- 
rotation band by a molecule at a given altitude. A single 
entry to Function XFER computes this number for radiation 
received fron molecules located either below or above this 
altitude. Equation (44) is summed over each altitude in- 
terval involved. 


Function XFER assumes that the values in RAD(1) and TTH(1) 
are the values corresponding to the lowest altitude over 
which the integration is to be performed. Thus, when 
called from Subroutine BANRAD to determine the effect of 
molecules above a given altitude with index #Z, the call- 
ing argument list has the for (VIBPOP(NZ,2), THR(NZ,2}, 
NALT, NZ, 2,KB(1)) for a band transitica between vibra- 
tional Levels 1 and 2 and band numbe. 1. 


Upon entry te Function XFER, the number of altitude inter- 
vals is found from the differance of NTOP and NBOT, and 
stored in NINT. If NTOP“NBOT, the value }.0 is returned 
by Function XFER. Otherwise, K is tested and control 
transferred to the appropriate part of fhe function. 


When K is one, AFER computes the radiative contribution 

to excitation from layers below the altitude level for 
which the computation is made. When K is two, the contri- 
bution from layers above is computed. 
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The computation fron a single level is found in each pas- 
Sage through the loop ending with Statement 120 (if K is 
one) or Statement 220 (if K is two). Before entry to 
either loap, RA is set to one, TB to zero, and EA3 to zero. 
RFAC, the ratio of ground state to upper state vibrational 
populations is computed and SUM is set to one-half. Within 
the appropriate loop, TB is incremented by the optical 
thickness obtained from array TTH. 1/2 Nyo(t) for the 
band is obtained from Function SLM and the value stored 
in EBS. The band function 1/2 L(t) ¥s computed from the 
value of Liz) obtained by a call to Function SL} and 

172 Nio(t) through use of the recurrence relation 


& L(t) = tft * Nyolt)] 


Equation (44) is than applied and the result added to the 
value ta SUM, 


HNENONIC VARIABLES 

EA3 temporary storage for value of the bard function by for 
the nearer of two adjacent integration points 

EB2 temporary storage for value of the band function 1/2 Nis 
for farther away of two adjacent integration paints 

£B3 tenpera:y storage for value of the band fanction Ly for 
the farther away of two adjacent integration points 

K integer parameter used to choose appropriate section of 
this function 

KB integer parameter used to choose appropriate sections of 
Functions SLML and Si 

Head indices in altitude array of top and botrom altitudes 


in layer of interes! 


Lod 


V1GPOP 


Non-Library 
Externals: 


falied by: 
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index used to select values from array VIBPOP 

(NTOP - NBOT), number of intervals 

(NINT#) ) 

index for altitude level 

index for altitude level 

relative vibrationa? population of never altitude point 
real array containing vibrational po;:lation of an upper 
state 

relative vibratioral population of ‘e-ther altitude point 
ratio of ground state population te apper vibrational 
Yevel population 

immediate result of integration 

optical distance from level for wrich calculation is being 
made to farther of two levels 

increment of optical thickness Latween two adjacent levels 
real array storing incrementa! «tical thicknesses 

rea) array storing vibrational »vepulations of ground and 
excited levels of a molecule 


SLMI, SL3 


BANRAD 
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Function SL3 


Purpose: 


z 
a 


Program Flow: 


Inputs: 


Output: 


This function yields the total abserption of radietion 

from & continuum source passing through a line absorbing 
medium. The result of this computation is used by Function 
XFER to compute the band function ue (t). 


SL3 computes the function aes) for parallel and per- 
pendicular transitions in linear molecules. The algorithm 
uses a polynomial approximation to values of (2) ccm 
puted for a reduced rotational temperature heB/kT = 0.002. 
The valules computed by SL3 are within five percent of 

the correct values for any reduced rotational temperature 
encountered in running the program. 


Upon entry, the absolute value of the argument TAU is de- 
fined as T. Some logic then is used to choose the appro- 
priate polynomial curve fit to the function which is 
computed and return is executed. 


K, TAU 


SLM1 


MNEMONIC VARIABLES: 


input parameter used to choose appropriate branch of the 
function, K equals 1 for parallel banc, K equals 2 for 
perpendicular band 

computed value of 1°(7) 

absolute value of argument of TAU 

input variable, optical depth ky 
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Non-Library 
Externals: 


| Called by: 
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argument used to compute curve fit 


natural logarithm of T 


None 


XFER 
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| 
| 


MNEMONIC VARIABLES: 


TAUMAX 


Non-i ibrary 
Externals: 


rotational constant (em) 

columnar number density of molecule (cm?) 
wavelength (ym) 

integer parameter 

integer parameter 

mass of molecule (gm/mole) 

band strength (on™  (atm-cm)~! 

the quantity hcB,/kT 

temperature (K) 

value of optical thickness returned by function 


SUMAX 


BANRAD 
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Function TAUMAX 


Purpose: 


Method: 


Program Flow: 


The purpose of this function is to compute the optical 
thickness of the strongest ‘ine in an infrared vibrational- 
rotational band. 


Function TAUMAX computes the optical thickness ko at the 
center of the strongest line of a vibrational-rotational 
band of a linear molecule, assuming Doppler line shape. 


The optical thickness at the center of a line with Doppler 
shape is 


a 38 1/4 Tine 
k, * 5.3612 x 10 Vr Nieto ty 


wnere M is the molecular mass (gm/mole), T is the abso- 
lute cmperarire (k), N le is the colum.ar number densiiy 
(molecules/cm zy 18 the line center wavelength (um) 
and Pind is the line strength in absorption (cm™? 
(atm-cm)~! at 30u x. Function TAUMAX computes the re- 
duced rotational temperature neB, /kT and transmits this 
as an argument to Function SJMAX which computes the frac- 
%1on of total hand streayth in the strongest line 

a isp TAUMAX then computes the optical thickness 


at the center of the strongest line, using ee in place 


of a in the above expression. 
KL, KB, RM, BV, T, FLAM, S, DEN 


TAUMAX 
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FLOWCHART: TAUMAX 


S1Ge] .438739*BV/T 


LM DERSSIAN(KL 8,516) | 


RETURN 


TAUMAX = (5.3613 x 10728{/RM og 
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Function SJMAX 


Purpose: 


Program Flow: 


Inputs: 


Qutput: 
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a RA RRS SW lle RY MSE ene Se eT soe ee cao eee Ae Se we cee Ste Se Pea 
Brace) aed Ba” a Qe eR ew ey 


The purpose of this function is to compute relative line 
strengths. 


Function SJMAX computes the relative strength of the 
strongest line in a parallel or perpendicular band of a 
linear molecule, assuming that all parallel bands are 
Ip aty transitions ami all perpendicular bands are 


1, os) transitions. The effect of negiecting the actual 
transition involved is negligitle for the uses of 
Function SJMAX. 


Functicn SJMAX computes the J values of the strongest 
line in the band, computes the strength from the appro- 
priate Hini-London formula, and divides this by the 
rotational partition function returned by Function QUPART. 
KL, KB, SIG 


SUMAX 


MNEMONIC VARIABLES: 


intermediate result in calculation 

parameter to chocse type of transition 

(KB = 1) -) transition) (KB = an - ) transition) 
paramete: passed to QJPART 

the quantity (heB, )/kT) 

strength of strongest line in band 

relotive strength of strovgest line in band 
effective rotational quantum number 
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Non-L ibrary 
Externals: 
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FLOWCHART: SJMAX 


START 


Lee 
nee 
120 120 


ran « H(Yarsie - 1 SUMAX = 0.0 


2 
A= 5 (210M + 1) A= TM +1 
RETURN 
@ 
200 
SUM = A*EXP(-SIG*TUM*(TUM + 1) 
-amticen 
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Function QUPART 


Purpose: The purpose of this function is to compute rotational par- 
tition fusctions for linear molecules. Refer to Subroutine 
BANRAD for definition of symbols. 


Method: Function QUPART computes the rotational partition function 
of a single vibrational level of a rigid linear molecule, 
assuming that the enerjy E) of a rotational state with 
quantum number J is 


Ee," heB J(u + 1) 


Provision is made for computing the rotctioral partition 
functions for vibrational levels in which rotational states 
with only even or odd J values are present. 


Program Flow: If K = 1, all rotational levels are assumed porulated. 
If K = 2, only even rotational levels are assumed popu- 
vated. uf K = 3, only odd rotational levels are assumed 
populated. 


The computation uses polynomial approximates to the 
partition function 


Q) = £)(20 + 1) exp (-E ,/kT) 
with SIG as the parameter. The range of applicability 
is limited to values of SIG less than 0.2, and the largest 
value of SIG encountered for the molecules presented in- 


cluded is 0.015. 


Inputs: K, °IG 
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Outputs: QUPART 


MNEMONIC VARIABLES: 


K parameter used to choose appropriate fit 

QUPART value of rotational partition function 

SIG the quantity (hoB /kT) 

X the value of SIG, defined for convenience in the curve 
fit 


Non-Library 
Externals: None 


Called by: SUMAX 


lll 


| 
| 


a3 


FLOWCHART: QUPART 


COMPUTE FROM A 
GIVEN CURVE FIT 


C) RETURN 
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USE A 
‘IFFERENT FIT 


RETURN 
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Rabati, 


THIRD FIT 


RETURN 


Function SIMP 


Purpose: 


Method: 


Program Flow: 


The purpose of this function is to compute an integral 
using Simpson's rule. 


Function SIMP computes the integral of a function between 
any two of three equally spaced points. Simpson's rule 
integration is used. 


Upon entry, the value of K is tested to choose the appro- 
priate branch of the function with the following results: 


K< 2. SIMP = (h/12) (8A (2) + SA (1) - A (3)) 
K= 2  SIMP = (h/12) (8A (2) + 5A (3) - A (1)) 
K>2 SIMP = (h/3) (4A (2) + A (1) + A (3)) 


Inputs: A, H, K 
Qutput: SIMP 
MNEMONIC VARIABLES: 
A array of values to be integrated 
H integration step size 
K parameter to choose method desired 
Non-Library 
Externals: None 
Called by: HANRAD, XAFER 


| tales « 
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FLOWCHART: = SIMP 


) START 


NEG POS 
K-2 


100 ZERO 


SImP=(@A(2) + SA(1) - A(3) Hy 


7) RETURN 110 


SIMP=(4A(2) + A(1) + A(3) epg 


() RETURN 120 


SIWP*LAA(2) + A(1) + AC3)Ie4 | 


RETURN 
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Non-Library 
Externals: 


Called by: 


argument used to compute curve fit 
natural logarithm of T 


None 
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Function SMi 


Purpose: 


Method: 


Program Flow: 


Inputs: 


Qutput: 


This function yields a factor needed to compute the true 
rate of absorption at an optical depth, ko» of radiation 
(from a spectrally constant source) in a resonant medium. 
Refer to Subroutine BANRAD for further explanation and 
definition of symbols used below. 


SM comoutes the function w(K.) for parallel and perpen- 
dicular transitions in linear molecules. The algorithm 
employs a polynomial approximation to values of "> (ko) 


computed at a reduced rotational temperature heB,/kT = 0.002. 


Upon entry the absolute value of the argument TAU is de- 
fined as T. Some logic then is used to chvose the appro- 
priate polynomial curve fit to the function which is com- 
puted and return is executed. 


K, TAU 


SM1 


MNEMONIC VARIABLES: 


SM1 


TAU 


input parameter used to choose appropriate branch of the 
function; K equals 1 for parallel band, K equals 2 for 
perpendicular band 

rate of absorption at optical thickness ko relative to 
zero optical thickness 


absolute value of argument Ti" 


<nput variable - optical depth ky 
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FLOWCHART: SM1 


SM1*POLYNOMIAL IN X 


TAU=ABS (TAU) 


OC RETURN 
YES 


Xs 5 (T-2) 
SM1=POLYNUMIAL IN X 


NO 


y 


110 ©) RETURN © RETURN oe 
P10 X=(T-7.0)/3.0 X=(T-7,0)/3.0 NO 
YES SM1=POLYNOMIAL IN X SMI=POLYNOMIAL IN X 
120 
© RETURN ©) RETURN YES 
220 
T>100 SMO} X2(9.21...)/2n(T)-3.0 X=(9,21...)/tn(T)-3.0 | NO 
SM1=(POLYNOMIAL IN X)/T SMl= (POLYNOMIAL IN X)/T 
— O RETURN © RETURN NES 
130] X=(9.21...)/2n(T)=-1.0 X=(9.21...)/£n(T)-1.0 1.230 
SMi=(POLYNOMIAL IN X)/T SM1=(POLYNOMIAL IN X)/T 


@ RETURN 


@ RETURN 
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Function SLM1 


Purpose: 


Progrem Flow: 


This function yields the factor used to obtain the actual 
absorption of radiation from an extended black vody sur- 
face passing through a resonant medium. f.efer to Subroutine 
BANRAD for further explanation and definition of symbols 
used below. 


SLM1 computes the function ae WP(k,) foy parallel and per- 
pendicular transitions in linear molecules. The algorithm 


uses a polynomial approximation to values of ge NP(k,) com- 
puted for a reduced rotational temperature heB, /kT = 0.002. 
The values computed by SLM1 are within five percent of the 
correct values for any reduced rotational] temverature en- 
countered in running the program. 


Upon entry the absolute value of the argument TAU is defined 
6s T. Some logic then is used to choose the appropriate 
polynomial curve fit to the function which {is computed and 
return is executed. 


K, TAU 


SLM] 


MNEMONIC VARIABLES: 


input parameter used to choose appropriate branch of the func- 
tion, K equals 1 for parallel band, K equals 2 for perpendicu- 
lar band 

computed value of ae NPR) 

absolute value of argument TAU 

input variable, optical depth ky 
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Non-Library 


Externals: 


Called by: 


argument used to compute curve fit 
natural logarithm of T 


None 


BANRAD 
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FLOWCHART : 


SLM 


TAUMABS ( 


X=(T-.5)/2 
SLM1=POLYNOMIAL IN X 


O RETURN 
110 
wo | X=(T-2.5)/1.5 
SLM1=POLYNOMIAL IN X 
YES RETURN 
120 
NO X=(T-7.0)/3.0 
T>10 SLM1=POLYNOMIAL IN X 
YES RETURN 
130 
NO a0 
X=(9.21...)/pT-3.0 
, SLMI= (POLYNOMIAL IN X)/T 
YES o 


RETURN 


X2(9.21...)/mnT-l. o) 
SLML= (POLYNOMIAL IN X)/ 


J RETURN 


X#(T-.§)/2 
SLM1=POLYNOMIAL IN X 


X=(T-2.5)/1.5 


X(9.21, 
)/t Saree 


START 


TAU 


RETURN 


SLM1*POLYNOMIAL IN X § 


SpetuRn 


X=(T-7.0)/3.0 
SLM1=POLYNOMIAL IN X 


RETURN 
 acerimasbtan ene eneemenenante iaetetatiandatietaneaateememmmamaiae Ueatmmbenen N 
X=(9.21,..)/ent-3.0 cas 
SLM1= (POLYNOMI) IN X)/T YT sai’ 
J aa YES 
po 
aera 1.0 
POLYNOMIAL IN x)/T 
RETURN 
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PROGRAM SPCTRA 


Purpose: 


Program Flow: 


. 


SPCTRA computes the spectral radiances in the 5 to 2 um 
interval corresponding to the band radiances computed by 
BCKGND. The input data to SPC7RA are the molecular band 
Tine constants used by, and vibrational population and 
band radiance data generated by SCKGND. In its present 
form, SPCTRA computes spectral radiances between 2.5 and 
25 um for limb views and endo-atmospheric zenith angles 
of 9 anc 90 degrees at 5 km intervals between 60 and 200 
km, 


When program execution begins, values are assigned to the 
integer variables NOUT, NLN and NALT and the Subroutine 
ROTATE, which reads in data for the molecular lines, is 
called. The outer DO range ending with Statement 500 is 
entered. The DO parameters are fixed at present to pro- 
duce two sets of spectra, using night and day input from 
Program BCKGND. The program may be made more flexible 

by reading in DO parameters prior to execution of the DO 
statement. The variable BLAM is set to 5.0 and the outrr 
DO range ending with Statement 105 is entered. The effect 
of the first two statements is to set up a table of wave 
numbers at 2 cm”! intervals between 2.5 and 25 um in real 
array AVLAM. These are the finest wave numbers at which 
spectral radiances may be printed out. The inner DO range 
ending with Statement 105 sets the values of output spec- 
tral radiances equal to zero. 


The DO range ending with Statement 110 sets temperature 
values in array TEMP to those corresponding to the U.S. 
Standard Atmosphere, 1962. Temperature values may also 
be read in prior to execution of the first DO statement. 
Following final execution of this DO range, temperatures 
are printed out. 
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Next, the number of molecular species is read in and stored 
in integer variable NSPTOT. The outer DO range ending with 
Statement 400 is entered and executed NSPTOT times. Within 
this DO range, molecular vibrational level and band data 
and vibrational population data are read in and mean temper- 
atures and total number of molecules per unit area in the 
vertical and horizontal directions are computed for each 
vibrational state. Total band radiances corresponding to 
limb viewing and vertical and horizontal endo-atmospheric 
viewing are read in. Relative line radiances are computed 
and normalized to the total band radiances. The spectra 
corresponding to the normalized line radiances are computed 
and stored in arrays SPECA, SPECB and SPECC. Details of 
these computations are given below. 


The code number of the species ISP is read in. The number 
uf vibrational levels NLEVEL, the number of bands NBAND and 
the molecular weight RM are read in. In the DO ranje ending 
with Statement 210, the level description codes LSC and LCC, 
the energy and the mean rotational constant for each vibra- 
tional level are read in. In the DO range ending with 
Statement 22. the band description codes LBC and LBU and the 
band strength for each band are read in. The number of alti- 
tudes NALT tor which vibrational populations are computed 
for the species and the number of altitudes NRLEV for which 
band radiances are computed by BCKGND are read in. Vibra- 
tional level populations are read in and stored in array CON. 


The mean temperatures and total number of molecules per square 
centimeter in the vertical and horizontal directions at 5 km 
intervals between 60 and 200 km are determined in the DO 

range ending with Statement 129. The product of temperature 
and vibrational level population is computed in the DO range 
ending with Statement 112. 
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The integrals of total number of molecules and the product of 
temperature and number density in the vertical direction are 
computed in the DO range ending with Statement 120. Simpson's 
rule integration is used. The total number of molecules in 

a given vibrational level above a given altitude is stored in 
array COLCNT. The corresponding mean temperatures, obtained 
by dividing the integral of the product of temperature and 
number density by the integral of number density, are stored 
in array AVTEMP, After exit from DO range, the column counts 
and temperatures are printed out. 


The DO range ending with Statement 128 performs a similar 
computation for colunn counts and mean temperatures in the 
horizontal direction. The integrations required are per- 
formed in Function QINT. Upon exit from this DO range, 
column counts stored in array HORCNT and mean temperatures 
stored in array BVTEMP are printed out. 


The inner 00 range ending with Statement 400 is now entered 
and executed for each band. This DO range controls the com- 
putation of the spectra. 


The altitude HT, the limb band radiance RADC, and endo- 
atmospheric radiances RADA (upward <:’rection) and RADS 
(horizontal direction) are read in the DO range ending with 
Statement 230. Band information is computed from band and 
level constants. A test is made on the species code ISP 

to determine whether the molecule is jinear or non-linear. 


Following the test of ISP, if the molecule is linear and 
is not in McClatchey's tabulation, line positions are com- 
puted in Subroutine ROTPOS. The DO range ending with 
Statement 250 is entered and executed for each of the 
twenty altitudes for which spectra are computed. Line 


1&5 


radiances are computed in Subroutine ROTSTR for vertical 
endo-atmospheric viewing. Successive calls to Subroutine 
FILTER compute the spectral contributions from lines in the 
P, Q and R branches. The sequence of calls to ROTSTR and 
FILTER is repeated for endo-atmospheric viewing in the 
horizontal direction and for limb viewing. Upon exit from 
this DO range, control is transferred to Statement 280. 


If the test on ISP determines that the molecular line data 
is available, contro] is transferred to Statement 250 and 
the DO range ending with Statement 270 is executed for each 
of the 28 altitudes for which spectra are computed. Spectra 
are generated from the line radiances by calls to Subroutine 
WATSTR and FILTER, for each of the three viewing conditions. 


Following the final execution of the DG ranges ending with 
Statement 400, the spectra are printed out in the order 
erdo-atmospheric viewing in the upward direction, endo- 
atmospheric viewing in the horizontal direction and limb 
viewing. 


Upon final execution of the DO range ending with Statement 
500, the program hi.lts with a STOP statement. 


Mnemonic Variables: 


AVLAM real array in which output wave numbers are stored 

AVTEMP real array in which mean temperatures in the vertical 
direction are stored 

BY real array in which rotationa! constants for each vibra- 
tional level are stored 

BVTEMP real array in which mean temperatures in the horizontal 


direction are stored 
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COLCNT 


CON 


RHOLAM 
HHOLSE 


HHOSTR 
HORCNT 


PRAD 


RADC 
RCM 
ROTEA 


real array in which total number of molecules per square 
centimeter above 60, 65, .... 200 km in each vibrational 
level are stored 

real array in which number densities for each vibrational 
level are stored 

real array in which wave numbers of band lines are stored 
real array in which lower state energies of band lines 

are stored 

real array in which strengths of lines are stored 

real array in which total number of molecules per square 
centimeter in each vibrational level in the horizontal 
direction at altitudes 60, 65, ... 200 kin are stored. 

real array in which line radiances are stored 

real array in which wavelengths of P branch of linear 
molecule are returned by Subroutine ROTPCS 

real array in which products of number density and 
temperature are stored 

real array in which P branch line radiances are returned 
by Subroutine ROTSTR 

real array in which Q branch line wavelengths are returned 
by Subroutine ROTPOS 

real array in which Q branch line radiances are returned 
by Subroutine ROTSTR 

real array in which R branch line wavelengths are returned 
by Subroutine ROTPOS 

real array into which vertical direction and endo-atmospheric 
band radiances are read 

real array into which horizontal direction endo-atinospheric 
band radiances are read 

real array in which limb view band radiances are read 

real array into which vibrational level energies are read 
real array used by Subroutines ROTATE and WATSTR, stores 
factors needed to compute line radiances in water vapor 
pure rotational spectrum 
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ROTLAM 


ROTLSE 


ROTST 


RS 


SPECA 


SPECB 


SPECC 


STR 


TEMP 


BLAM 


CENTER 
HT 


H1, H2, H3, 
H4, H5 


RBASE 


real array in which water vapor pure rotational line 
wavelengths are stored 

real array in which lower state energies of water vapor 
pure rotational lines are stored 

real array in which absorpticn strencths of water vapor 
yure rotationa’ lines are stored 

real array in which R branch line radiances are returned 
by Subroutine KOTSTR 

real array in which spectral radiances for vertical view- 
ing endo-atmospheric geometry are generated 

real array in which spectral radiances for endo-atmospheric 
horizontal viewing are generated 

real array in which spectral radiances for limb viewing 
are generated 

real array used to store band strengths 

real array used to store atmospheric temperatures 

real variable used as temporary storage during generation 
of values in array AVLAM 

real variable used to store wave number of band center 
real variable into which altitudes at which band radi- 
ances are computed by BCKGND are read 

real variables used to store distance increments used 

in computing total column counts and mean temperatures 

in the horizontal direction 

real variable used to store distance of base altitude 
from center of earth in computing horizontal column 
counts 

real variable used to store value of RBASE**2 

real variable used tu store value of molecular weight 
real variable used to store distance from center of 
earth of given altitude leve) 

real variable used to store current value of band strength 
real variable used as temporary storage in computing 
horizontal column counts 
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SVH real variable used as temporary storage in computing hori- 
zontal column counts 


TC real variable used to accumulate integral of number density 
in computing vertical column counts 

TCH real variable used to accumulate integral of number density 
in computing horizontal column counts 

TP real variable used to accumulate integral of product of 
temperature in the vertical direction 

TPH real variable used to accumulate integral of product of 


temperature and number density in computing mean tempera- 
ture in the horizontal direction 


LBC integer array scoring level numbers of initial and final 
states of a band 

LBU integer array storing band type and symmetry information 

LCC integer array storing vibrational level information 

LSC integer array storing vibrational level information 

I integer variable used as index in DC loops 

ISP integer variable into which species identification code 


is read (2 for CO,, 3 for H 
6 or 7 for ozone) 

J integer variable used as index in DO loops 
inteyer variable (set equal to 28 in current program) 
used as final parameter in DO statement controlling com- 
putation of column counts and mean temperatures. The 
value of K is the number of altitudes for which spectra 


20> 4 for NO, 5 for No0 and 


are computed 

KB integer variable used to store band type code 

KK integer variable used as final parameter in DO range end- 
ing with Statement 125, computed from K and I 

Kt integer variable used to store band symmetry code 

NA integer variable used as altitude level index in com- 
puting horizontal column counts 

NALT integer variable storing total number of altitude levels 
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NB 


NBAND 


NBC 


NCASES 


NF 


NI 


NL 


NLEVEL 


NLN 


NRLEV 


NSP 


NSPTOT 


NZ 
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integer variable used as index for number of current 

band 

integer variable storing total number of bands considered 
for a given species 

integer variable used to store current value of LBC 
integer variable used as index for spectrum set 

integer variable used as index for lower vibrational state 
of band 

integer variable used as index for upper vibrational state 
of band 

integer variable used as index for vibrational level 
number 

integer variable in which is stored total number of vibra- 
tional levels for a molecular species 

integer variable controlling number of lines treated in 
each branch of a band of a linear molecule, set to 100 

in current program 

integer variable into which is read number of altitudes 

at which band radiances for a given species are computed 
by BCKGND 

integer variable used as index for DO range executed 
NSPTOT times 

integer variable into which is read total number of species 
to be included in a spectrum 

integer variable used as altitude index in DO loops 
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FLOWCHART : SPCTRA 


READ SPECIES CODE NUMBER | 


DO LOOP 210 READS VIBRATIONAL LEVEL DATA 
DO LOOP 220 READS BAND CONSTANTS 


DO LOOP 221 READS VIBRATIONAL POPULATIONS 


SET NL=1 FOR FIRST LEVEL 


DO LOOP 112 FORMS PRODUCT OF TEMPERATURES 
AND VIBRATIONAL POPULATION 
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DO LOOP 120 FORMS INTEGRAL OF PRODUCTS FOR VERTICAL DIRECTION 


DO LOOP 128 FORMS INTEGRAL OF PRODUCTS FOR HORIZONTAL DIRECTION 


: © 


YES 


IF SPECIES IS WATER, CALL ROTATE, WATSTR 
AND FILTER TO 00 CALCULATIONS FOR ROTATIONAL LINES 


SET NB=1 FOR FIRST BAND 


DO LOOP 230 READS BAND RADIANCES FOR ALL LEVELS 


CALL ROTATE, WATSTR AND FILTER TO PERFORM SPECTRAL CALCULATIONS 


ene 
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MAKE CALLS TO LINOUT TO PRINT SPECTRAL VALUES 
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Subroutine FILTER 


Purpose: FILTER computes the contribution to spectral radiance trom 
2ach line in a band or sub-band, in the interval 2.5 um 
to 25 um. Provision is made for varying resolution in 
either wavelength or wave number space. Sufficient stor- 
age is provided to compute radiances at two wave number 
intervals with triangular slit resolution of five wave 
numbers. 


Method: FILTER searches an array of emission Vine strengths and 
assigns a weighted strensth to an output wave number 
or wavelength. 


Program Flow: The outer DO range ending with Statement 150 controls the 
distribution of spectral radiance from the Ith of N given 
lines, and computes the indexes JS and JF of the first 
and last output points in array SPEC to which a given 
line contributes. The inner DO range ending with State- 
ment 120 computes the contribution of a given line to 
spectral radiance at output wavelengths or wave numbers. 


Inputs: RAD, WVLTH, N, KOMPUT 
Output: SPEC 


MNEMONIC VARIABLES 


SPEC real array into shich spectra radiances are added at reg- 
ular intervals between 2.5 and 25 um or 400 and 4000 
wave numbers. 
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RAD real array containing radiances of individual lines 
WULTH real array containing wavelengths of individual! lines 
194 


¢ peer ny ETN 
Sain ; 


KOMPUT 


RADT 


RLAM 


uF 


JS 


o¢ 


integer variable containing number of lines in band or 
sub-band for which spectral radiance is to be computed 
integer variable which determines number of output spectral 
points and whether wavelength or wave number output is 
made 

real variable used as temporary storage for radiance 
value of a single line 

real variable used as tempurary storage for wavelength 
value at which spectral radiance ts computed 

real variable used as temporary storace for wavelength 
of a single line 

integer variable used as index for outer D0 Joop which 
treats each line in turn 

interger variable used as index for inner DO loop which 
assigns spectral contribution of each line to output 
wavelength 

integer variable used as final parameter in inner 00 
loop, index of greatest output wavelength for a given 
line 


integer variable used as initial parameter in inner 00 
loop, index for shortest output wavelength for a given 
line 
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Purpose: 


Program | low: 


Subroutine ROTATE 


Subroutine ROTATE is used to read in spectral line data 
from the compilation of McClatchey, et al. The {nput 
data consists of Tine positions and strengths for vibra- 
tional - rotational bands of water, carbon dioxide, ozone 
and nitrous oxide, and for the pure rotational spectrum 
of water, This data is processed for use by other 


subroutines. 


The DO range ending with Statement 110 reads in data on 
the strongest lines in a vibration - rotation band. Tints 
data consists of the wave number of the Tine, the Tine 
strength at 296 K and the wave number of the lower rota- 
tianal state. ROTATE converts the wave number to wave- 
Teneth in pm stored in array HHOLAM, converts the Tine 
strength to units cw atmo dt P73.35 & stored in array 
HHOSTR and stores the wave number of the Jower state in 


array HHOLSE. 


The rotational partition function for water vapor is com 
puted next if water vapor data is being read. the 0 
range ending at Statement 230 ts used only tor water vapor 
pure votationa? Lines and computes the waveiength (stored 
in array ROTLAM), Finstein A coeffieiept (stored fa array 
tA), relative line strenuth in absorption at Ofa.i & 
(storad tn array ROTST), relative Tine strenath on etins 
sion at O73.E8 & (stored dn array ROTEAY and the wave 
nuber of the lower rotational state (stored mn arrey 
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MNEMONIC VARIABLES 


EA 


EAL 


EAL" 


ENERGY 


FNL 


ENV 


HHOLAM 


HHOLSE 


HHUSTR 


ROTEA 


ROTLAM 


ROTLSE 


ROTST 


WVLTH 


LAFAC 


QR 


real variable used to store Einstein A coefficient for 
water vapor pure rotational spectrum line 

real variable used to store intermediate results in 
computing pure rotational line stvengths in emissic. 
real variable used to store intermediate results is com- 
puting pure rotational line strengths in absorptio: 

real variable into which is read wave numbers of sire 
rotational lines 

real variable into which is read wave numbers of ‘ower 
states of pure rotational lines 

real variable into which is read wave number of .yper 
state of pure rotational line 

real array in which wavelengths of vibration rotation 
band ines are stored 

real array in which wave numbers of lower st«tes of vibra- 
tion - rotation band are stored 

real array in which line strengths of vibration - rotation 
band are stored (om e-pom”! STP) 

rea) array in which pure rotational line strengths in 
emission (273.15 k) are stored 

rea! array in which pure rotational tine wavelengths are 
stored 

real array in which lower state wave number of pure 
rotational lines are stored 

real array in which pure rotational line strengths in 
absorption are stored 

real array in which wavelengths of vibration - rotation 
band lines are stored (redundant) 

“eal variable used to store pre-factor in computing 
finstein A coefficients for pure rutational lines 

real variable in which is stored approximate rotational 
partition function for water vapor molecule 
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QRFAC real variable into which is read pre-factor for computing 
7 rotational partition function 
7 RA real variable into which is stored first water vapor 
rotational constant 
RB reai variable into which is stored second water vapor 
rotational constant 
RC real variable into which is stored third water vapor 
rotational constant 
TEFAC real variable in which is stored exponential factor used 
in converting line intensities at 296 K to line strengths 
at 273.15 k 
TOFAC real variable in which is stored factor reiating rotational 
function at 296 K to that at 273.15 k 
I integer variable used as index in DO loop 
NLINES integer variable into which is read number of lines in- 


cluded in pure rotetional spectrum and in any single 
vibration - rotation band 
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Subroutine WATSTR 


Purpose: 


Program Flow: 


WATSTR computes the radiances of individual lines in vibra- 
tion - rotation bands and in the pure rotational spectrum 
of water vapor. The radiances are computed with the 
assumption that the atmosphere in the line-of-sight may 

be treated as having uniform temperature. The radiance, 
computed for the vibration - rotation bands are scaled 

to match the total band radiances produced by the BACK- 
GROUND Program. Absolute radiance values are computed 

for the water vapor pure rotational spectrum. 


Upon entry to WATSTR, a conversion factor TQFAC is com- 
puted for later use in computing the effect of charging 
rotational partition function with temperature. Control 
is then transferred to the appropriate part of the sub- 
routine as determined by the value of N. 


For computation of the radiances of lines in the water 
vapor pure rotational spectrum (N=1), control is trans- 
ferred to Statement 100. A radiance factor RADFAC and 
an optical thickness factor TAUFAC are computed. Then 
the optical thickness and radiance values corrected for 
optical thickness and induced emission are computef for 
each line in the DO range ending with Statement 110. 
Control is returned to the main program, Line constants 
generated by Subroutine ROTATE are available in the arrays 
of common area ROTCOM. Radiance values are returned in 
array HRAD. 


Control is transferred to Statement 200 fur computation 
of radiances of lines in the other bands. An optical 
thickness factor TAUFAC and factor TEFAC used in comput- 
inu populations of rotational states at temperature T 
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HH()!.AM 
HHOLSE 
HHOSTR 
HRAD 


ROTEA 
ROTLAM 
ROTLSE 
ROTST 
HSTR 


are computed. The variable SUM, used to accumulate the 
sum of relative radiances, is set to zero. In the Dd 
range ending with Statement 210, the otpical thickness 

of each line is computec and the relative radiance value, 
corrected for cptical thickness, is computed and stcred 
in array HRAD. The sum relative radiances is accumulated 
in SUM. Upon exit from this DO range, the total radiance 
RAD is divided by SUM and stored in RADFAC. In the DO 
range ending with Statement 220, absolute radiance values 
are obtained and stored in array HRAD. Control is then 
returned to the main program. 


MNEMONIM VARIABLES 


real variable containing total band strength for vibration 
- rotation band 

real variable containing molecular weight of molecule 
real variable containing value of temperature 

real variable containing total number of ground state 
molecules per square centimeter in line-of-sight 

integer variable determining whether entry to WATSTR re- 
sults in computation for water vépor pure rotaticnal 
spectrum (N=1) or other band (N=2). 

real array, values as in ROTATE 

real array, values as in ROTATE 

real array, values as in ROTATE 

real array used to return radiance values of individual 
lines 

real array, values as in ROTATE 

real array, values as in ROTATE 

real array, values as in ROTATE 

real array, values as in ROTATE 

real variable used as temporary storage in computing 
optical thickness 


one ee 
¥ 


ae Seen aie Sarerennn anh re ee 


Mercere Fe ns ce rire 


PEATE PLOT ITT BAS IESG EOE RTE PT OTT, BP PRT LP PA SE TT ET OT TE CE Es OE Oe, la IB a a ig Lae Pe 


ADFAC 


SUM 


TAU 


TAUFAC 


TEFAC 


TQFAC 


rea} variable tn which Js stored pre-factor used in com- 
puting radiance values 

rea] variable in which partial sums uf relative radiances 
are accumulated 

real variable in which optical thickness at line center 
is generated 

real variable in which pre-factor fer optical thickness 
is stored 

real variable in which factor used tn computing rotational 
populations is stored 

real variable in which factor converting rotational par- 
tition function js stored 
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Subroutine ROTSTR 


Purpose: 


Program Flow: 


Subroutine ROTSTR computes the relative radiances of lines 
in the vibration - rotation bands of those linear molecules, 
which are not included in the line parameter compilation 

of McClatchey, et al. The radiances are corrected for ef- 
fects uf optical thickness, and absolute values then de- 
termined by normalizing to total) band radiances computed 

by the BACKGROUND Program. As for non-linear molecules, 

a single temperatures 13 assumed for tne molecules in the 
line-of-sight. 


Upon entry to ROTSTR, the reduced rotational temperature 
SIG is computed, relative strengths of the P. Q and R 
lines for J=0 are assigned and values assigned to the in- 
crement of J values DJ and tre variable J value AJ. De- 
pending on the value of NM, relative strengths of rota- 
tional lines are computed in ether the DO range ending 
with Statement 110 (paralte! transition) or the DO ranae 
ending with Statement 120 (perpendicular transition). 
Control is transferred to Statement 200 and the pre-factor 
TAUFAC used in computing optical thickness °s evaluated, 
The variable SUM which accumilates the sum or relative 
radiances is set to zero. In the 00 range ending with 
Statement 310, optical thicknesses for lines in P, Q and 
R branches and retutive radiances for these lines, cor- 
rected for the effects of optical thickness, are computed. 
In the 00 range ending with Statement 320, absolute line 
radiances normalizea to the total band radtance RAD are 
computed, Control is then returned to the moti program. 
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MNEMONIC VARIABLES 


AJ 


DJ 


RS 


NF 


real variable containing total band strength 

real variable containing molecular weight of molecule 
real variable containing mean rotational constant of 
lower vibrational state 

real variable containing varying value of rotational 
quantum number J of lower vibrational state 

“eal variable containing value of increment of rotational 
quantum number 

real variable containing mean temperature of emittina 
laver 

real variable containing number of molecules in lower 
state per square centimeter in line-of-sight 

real variable containing band strength at 300 K 

real array in which radiances of lines in P branch ae 
revurned 

real array in which radiances of lines in Q branch are 
returned 

real array in which radiances of Jines in R branch are 
returned 

integer variable determining whether pand is parallel 
(NP-1) or perpendicular )NP=2) 

mteger variable determining whether all lines are pres- 
ent (N01) or alternate lines are missing (NO-2} 

real array containing wavelengths of lines in P branch 
real array containing wavelengths of Jines in @ branch 
real array containing wavelengths of lines in R branch 
res] variable used to store exponential factor ia computa- 
tion of Tine strengths 


real variable storing optical thickness in a P branch line 
real variable storing optical thickness ir a 0 branch line 
real variable storing optical thickness of R branch ine 


SIG 
SUM 


TAUFAC 


real variable storing reduced rotational temperature 

real variable which accumulates sum of relative line 
radiances 

real variable storing pre-factor used in computing optical 
thickness at line centers 
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Subroutine ROTPOS 


Purpose: 


Method: 


Program Flow: 


ROTPOS computes the wavelength of tines in the P, Q and 

R branches of the infrared vibrational - rotational bands 
of linear molecules neglecting spin and A or ¢ doubling. 
It is used for bands not included in the compilation of 
McClatchey, et a}. 


The line positions are computed using Equations (IV, 19), 
(1V, 20) and (IV, 22) of Herzberg (1945). It should be 
noted that use of these simple formulae results in an 
error in line positions which may approach 0.1 ym for 
large J values in transitions between doublet states. 


Upon entry to ROTPOS, the auxiliary values BDIF, 8SUM, 
RCENT and TRI used in the computation are computed from 
the band center position and the rotational constants. 
The initial J value is set to zero and provisions made 
tor omitting alternate lines by setting DJ equal to NO. 
This latter starting condition means that the alternate 
lines omitted are always those with odd J values. The 
wavelengths of 100 lines in each branch in units po are 
then computed in the DO range ending with Statement 110, 
The wavelengths are printed out and control returned to 
the main program, 


MNEMONIC. VARIABLES 


AJ 


ROLF 


8SUM 


real variable in which current rotational quantum number 
J is stored 
real variable storing difference of rotational constants 
RVP and BVPP 


real variable storing sum of rotational constants 


AL) 


DJ 


RCENT 


TRI 


CENTER 
BYP 


BVPP 


ND 


real variable storing increment of rotational] qu. ‘tum 
number 

real variable containing constant used in computing R 
hranch line positions 

real variable containing constant used in computing R 
branch line positions 

real variable giving the wave number of the band canter 
real variable giving the mean rotational constant of the 
upper vibrational state 

real variable giving the mecn rotational constant of the 
lower vibrational state 

real array in which wavelengths of the P branch are re- 
turned 

real array in which wavelengths of the Q branch are 
returned 

real array in which wavelengths of the R branch are 
returned 

integer variable determining whether all lines ina 

band are present (ND=1) or whether alternate lines ave 
missing (ND=2) 
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